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SUMMARY 

The  soils  used  in  this  study  were  a  sand  and  two  sand- clay 
mixtures.   Kneading  compaction  was  used  to  prepare  samples  of  sand- 
clay  mixtures  while  pluvial  compaction  was  used  to  prepare  compacted 
samples  of  sand.   The  testing  program  consisted  of  the  permeability 
test  and  the  pore  size  distribution  (PSD)  test.   A  technique  for  the 
preparation  of  PSD  specimens  of  cohesionless  sands  has  been  developed. 
The  mercury  intrusion  technique  was  used  to  conduct  PSD  test.   Falling 
head  tests  under  back  pressure  were  used  to  measure  permeability. 

Mathematical  interpretations  of  PSD  were  given  in  detail.   For 
the  soils  studied,  both  the  pore  size  distributions  and  pore  size 
density  functions  were  presented.   Pore  size  density  function  of  sand 
studied  showed  a  single  modal  characteristic  on  a  log  diameter  scale 
while  that  of  sand-clay  mixtures  showed   bimodal  characteristics 
on  a  log  diameter  scale.   The  influence  of  varying  water  content  and 
compactive  effort  on  the  fabric  of  the  sandy  soils  studied  was  charac- 
terized by  their  pore  size  distributions  and  pore  size  density  func- 
tions . 

A  permeability  model  which  related  the  pore  size  distribution  to 

the  permeability  of  soils  has  been  proposed.   The  prediction  of  per- 

2 
meability  by  the  proposed  model  was  excellent  (e.g.  R  =  0.98)  for 

3 
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both  the  author's  and  Garcia-Bengochea's  (1978)  data.  In  other  words, 
the  proposed  permeability  model  can  provide  an  excellent  predictive 
tool.   Thus  a  functional  relationship  between  the  PSD  and  the  permea- 
bility of  compacted  soils  is  clearly  demonstrated. 


CHAPTER  I   INTRODUCTION 

1.1   FABRIC  AND  ITS  ROLE  IN  GEOTECHNICAL  ENGINEERING 

The  characterization  of  soil  fabric  is  a  very  important  aspect  of 
geotechnical  engineering.   It  can  lead  to  better  understanding  and  pre- 
diction of  engineering  properties  of  soils.   The  term  "soil  fabric" 
has  widely  been  used  without  a  common  definition,  e.g.,  Terzaghi  and 
Peck  (1948),  Mitchell  (1956),  Scott  (1963),  and  Oda  (1972).   However, 
the  following  definition  given  by  Brewer  (196A)  is  preferred: 

"Soil  fabric  is  the  physical  constitution  of  a  soil  material 

as  expressed  by  the  special  arrangement  of  the  solid  particles 

and  associated  voids." 

Studies  at  Purdue  University,  e.g..  Diamond  (1970,  1971),  Bhasin 
(1975),  Reed  (1977)  and  Garcia-Bengochea  (1978)  among  others,  have 
shown  that  the  fabric  of  compacted  clayey  soils  can  be  identified  by 
their  pore  size  distribution  (PSD)  and  that  some  engineering  properties 
of  clayey  soils  can  be  quantitatively  correlated  with  their  PSD.   More- 
over, the  engineering  properties  of  sandy  soils  are  also  believed  to 
be  affected  by  their  fabric  characteristics,  e.g.,  Oda  (1972,  1976), 
Mahmood  (1973),  and  Mitchell,  et  al.  (1976).   Figure  1-1  illustrates 
the  concept  of   sand  fabric.   Two  different  arrangements   of  particles 
could  have  identical  void  ratios  and   grain  size   distributions, 
but  their  compressibility  characteristics  would  be  obviouly 
different  (Leonards,  as  cited  by  Holtz  and  Kovacs,  1981).   Similarly, 
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the  permeability,  liquefaction  properties,  and  the  static  strength  of 
sands  are  strongly  affected  by  sand  fabric. 

1.2   OBJECTIVE  OF  THIS  STUDY 

This  study  is  part  of  a  research  project  entitled  "The  Effects  of 
Pore  Size  Distribution  on  Permeability  and  Frost  Susceptibility  of 
Selected  Subgrade  Materials".   The  specific  objectives  are: 

1)  to  determine  the  fabric  characteristics  of  compacted  sand 
and  sand-clay  samples  by  their  PSD  as  measured  by  the  mercury 
intrusion  technique;  and 

2)  to  develop  a  theoretical  permeability  model  which  correlates 
the  permeability  with  the  PSD  of  the  selected  soils. 


CHAPTER  II   THEORETICAL  BACKGROUND 

2.1   PORE  SIZE  DISTRIBUTION 

2.1.1  Theoretical  Background  of  Pore  Size  Distribution  and  Pore 
Size  Density  Function 

The  pore  system  has  an  extremely  complex  geometry.   Therefore,  to 

characterize  a  soil  fabric  by  means  of  pore  size  and  its  distribution 

is  a  rather  abstract  concept.   In  fact,  there  is  no  clear  definition 

of  pore  size  and  therefore  its  distribution.   In  the  mercury  intrusion 

technique  (see  Section  2.1.2)  a  cylindrical  pore  shape  model  was  assumed, 

i.e.,  the  circular  diameter  was  considered  to  be  its  corresponding 
pore  size. 

Based  on  Washburn's  (1921)  equation: 

P  =  -  ^Tc£s0  (2-1) 

d 

where  P  =  absolute  pressure  being  applied 
T  =  surface  tension  of  mercury 

0  =  contact  angle  between  mercury  and  the  pore  wall 
d  =  the  size  of  pores  which  can  be  intruded  at  the 
pressure  P;  called  the  apparent  pore  diameter. 
With  the  assumption  of  constancy  of  the  surface  tension  and  contact 
angle.  Equation  (2-1)  implies  that  as  pressure  increases,  mercury  will 
be  intruded  into  progressively  smaller  pores.   For  each  pressure  incre- 
ment, mercury  is  forced  into  the  accessible  soil  pores  of  a  diameter 


larger  than  or  equal  to  that  which  is  calculated  by  Equation  (2-1)  at 
the  corresponding  pressure.   Meanwhile  the  cumulative  volume  of  pore 
space  at  any  corresponding  pressure  is  also  recorded.   Then  a  PSD  can 
be  generated  (see  Equation  2-6d  and  its  interpretation). 

Since  there  is  no  clear  physical  meaning  of  pore  size,  Scheidegger 
(1957)  proposed  a  probabilistic  approach  to  define  a  pore  system.   With 
slight  modification  to  Sche idegger 's  ingenious  idea  and  the  pioneering 
work  of  Drake  and  Ritter  (1945)  and  Garcia-Bengochea  (1978) ,  the  PSD 
and  pore  size  density  function  can  be  derived  as  follow  : 

Assume  X  is  a  random  variable  denoting  a  pore  diameter  at  any 
point  within  the  pore  space.   It  is  the  diameter  of  the  largest  sphere 
which  contains  this  point  and  remains  wholly  within  the  pore  space. 
Thus,  a  "pore  diameter"  can  be  assigned  to  each  point  in  the  pore  space, 
and  a  PSD  can  be  defined  by  determining  what  fraction  of  pore  space  has 
a  pore  diameter  between  x  and  x  +  dx.   ( Scheidegger ,1957) 

Pore  size  distribution  can  be  expressed  as 


F  (x)  =     f  (x)  dx  (2-2a) 

J-x. 

where  F  =  PSD  of  the  soil  being  considered 

f  =  the  probability  density  function  of  pore  size 

and 


/    f^(x)  dx  =  1  (2-2b) 

Jo 

The  probability  of  pores  with  diameters  ranging  from  x  to  x  +  dx 

is  f(x)  dx.    Now  denoting   the  volume  of  pores  having  diameters  be- 
tween X  and  X  +  dx  as  dV,  we  have 


dV  =  s(x)  •  f  (x)  dx  (2-3) 


where  s(x)  =  a  shape  factor  and  is  a  function  of  x  (i.e.,  it 

depends  on  the  pore  shape  model  chosen) . 
Recalling  Equation  (2-1)  and  noting  that  the  diameter  is  x,  we  have 

Px  =  constant  (2-Aa) 

or 


dx  =  -  ^  dP  (2-4b) 


Substituting  into  Equation  (2-3)  ,  we  have 


dV  =  -s(x)  •  f  (x)  .  dP/P  .  X  (2-5a) 


or 

P       d(V   -  V) 

f  (x)  =  ^  •  ~- (2-5b) 

X      X  •  s(x)       dP 

where   V  =  total  void  volume  that  can  be  intruded  by  mercury 

V  -  V  =  a  measurement  of  pore  volume  with  diameter  equal  to  or 
greater  than  x 
P  =  corresponding  pressure  needed  to  intrude  mercury  into 
pores  with  diameter  x. 
The  terms  x,  f  (x) ,  s(x)  are  the  same  as  defined  previously  in  Equations 
(2-2)  and  (2-3) .   Note  that  the  total  void  volume  that  can  be  intruded 
by  mercury  is  not  necessarily  the  same  as  that  existing  in  the  specimen. 
Therefore  the  pore  size  density  functions  based  on  different  defini- 
tions of  V  are  different  from  each  other  (see  Section  4.3.2).   Also 
note  that  s(x)  is  the  problem  in  determining  f  (x)  of  Equation  (2-5b)  . 
One  way  to  treat  it  is  to  introduce  the  volume  pore  size  density  func- 
tion, f   (x)  as  follows: 
vx 


f   (x)  =  s(x)  •  f  (x)  (2-6a) 

VX  X 

Substituting  Equation  (2-5b)  into  Equation  (2-6a)  ,  we  have 

p   d(V  -  V) 

f   (x)  =  -  .  ^^r^ (2-6b) 

VX      X       dP 

Substituting  Equation  (2-4b)  into  Equation  (2-6b)  ,  we  have 

d(V^  -  V) 

f   (x)  =  -  ^ (2-6c) 

VX  dx 

Integrating  Equation  (2-6c)  gives 

F   (x)  =  V^  -  V  (2-6d) 

VX         t 

where  F   (x)  =  the  cumulative  intruded  volume  which  is  directly  re- 
vx 

corded  in  the  mercury  intrusion  technique. 

So  far,  the  above  derivations  and  definitions  provide  the  theoret- 
ical development  of  pore  size  distribution  and  its  density  function. 
For  any  pore  diameter  which  can  be  calculated  from  Equation  (2-1) 
with  each  corresponding  pressure  being  recorded,  F   (x)  can  be  ob- 
tained from  the  volume  intrusion  record.   Hence  a  PSD  Pl.e.,  F   (^) 
vSc  x3  can  be  generated.   The  pore  size  density  function,  f   (x) ,  can 
then  be  determined  from  Equation  (2-6c)  . 

The  pore  size  distribution  defined  above  is  in  terms  of  pore  volume. 
In  many  cases,  it  is  convenient  to  express  the  PSD  in  terras  of  pore 
volume  per  gram  of  specimen.   It  is  also  a  standard  presentation  of  the 
PSD  test  as  specified  by  an  ASTM  standard  being  developed  by  Subcommittee 
D18.14.   Moreover,  some  researchers  prefer  to  express  the  PSD  in  dimen- 
sionless  terms  (e.g.,  a  percentage)  for  the  purpose  of  convenient  math- 
ematical treatment.   The  following  expressions  provide  a  theoretical 


basis  for  this  consideration. 

p      d(V  -  V) 

where  f   (x)  =  density  function  of  pore  diameter  with  percentage  as 

its  ordinate 
and 


F   (x)  =  /    f   (x)  dx  (2-7b) 

px     J^  px 

where  F   (x)  =  cumulative  distribution  of  pore  size  (for 
px 

diameter  >^  x)  with  percentage  as  its  ordinate; 
this  is  simply  called  percent  greater  by  volume. 
Furthermore,  Equation  (2-6b)  can  also  be  modified  to  the  follow- 
ing expression: 

p  '     d(V   -  V) 
s 

where  V  =  volume  of  soil  specimen 
s 

f   (x)  =  density  function  of  pore  diameter  with  a  relative  frequency 
scale  such  that 


/    f   U)  dx  =  n 
/    nx 

where  n  =  porosity  of  specimen. 

The  multiplication  of  f   (x)  by  a  constant  c  (see  Section  4.3.2) 

nx     -^ 

is  called  "porosity  frequency"  by  Garcia-Bengochea  (1978)  and  White 
(1980). 

A  detailed  discussion  of  the  presentation  of  the  PSD  test  results 
on  the  soils  studied  is  given  in  Section  4.3.2. 


2.1.2   The  Measurement  of  Pore  Size  Distribution  and  Its  Limitation 
The  measurement  of  the  PSD  of  a  soil  is  a  technique  which  allows 

the  engineer  to  indirectly  examine  the  soil  fabric.   Methods  available 

for  the  determination  of  PSD  are  as  follows: 

a)  Indirect  determination  from  representative  grain  size  distribu- 
tions by  applying  probabilistic  theory  (e.g.  Muruta  and  Sato, 
1969); 

b)  Indirect  determination  by  scanning  electron  microscopy  (SEM) 
techniques  with  probabilistic  theory; 

c)  Capillary  suction  technique; 

d)  Mercury  intrusion  technique. 

Detailed  discussion  of  all  these  techniques  except  (a)  are  given  by 
Bhasin  (1975)  and  Lowell  (1979).   Bhasin  (1975)  had  also  pointed  out 
that  the  mercury  intrusion  technique  is  the  most  suitable  method  for 
the  measurement  of  the  PSD  of  soils. 

The  concept  of  forcing  mercury  into  material  pores  to  determine  their 
sizes  was  first  proposed  by  Washburn  (1921) .   This  led  to  the  develop- 
ment of  mercury  intrusion  technology  which  is  based  on  the  principle 
that  the  surface  tension  of  a  non-wetting  liquid  will  oppose  the  entry 
to  the  liquid  into  a  small  pore  of  a  material.  Washburn  stated  that  this 
opposition  could  be  overcome  by  external  pressure,  and  that  the  exter- 
nal pressure  required  was  inversely  proportional  to  the  pore  diameter, 
as  indicated  by  Equation  (2-1) . 

Perhaps  the  primary  limitation  of  the  mercury  intrusion  technique 
is  that  it  does  not  always  measure  the  diameter  of  the  actual  pore 
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but  rather  the  diameter  of  the  channel  leading  into  the  pore,  which 
obviously  could  be  smaller.   As  a  result  of  "ink-bottle-pore  effect", 
discussed  by  Ritter  and  Drake  (1945),  Orr  (1970),  and  Ahmed  (1971) 
and  the  other  Purdue  researchers,  the  measured  pore  size  is  called 
the  limiting  pore  diameter  or  called  apparent  pore  diameter  (e.g.,  by  a 
proposed  ASTM  standard  D18.14).   An  attempt  has  been  made  by  Cebeci, 
et  al.  (1978)  to  evaluate  the  arrangement  of  ink-bottle  pores.   How- 
ever, there  is  still  no  way  to  obtain  the  "true"  pore  size  distribution. 
This  problem  is  of  minor  impor''ance  in  the  case  of  flow  through  porous 
materials,  since  the  pore  channels  the  water  has  to  pass  through  are 
probably  equivalent  to  those  for  mercury. 

The  second  problem  with  mercury  intrusion  technology  is  that  cy- 
lindrical pores  are  assumed.   Bhasin  (1975)  pointed  out  that  this 
■assumption  may  be  fairly  reasonable  since  the  difference  in  calcu- 
lated PSD's  due  to  different  assumptions  of  pore  shapes  are  relative- 
ly small  (i.e., within  an  order  of  magnitude)  when  compared  with  the 
range  of  pore  sizes  that  are  possible  (i.e., up  to  five  orders  of  mag- 
nitude) . 

Another  limitation  is  that  mercury  intrusion  porosimetry,  in 
common  with  many  other  methods,  is  only  capable  of  intruding  pores 
that  are  open  to  the  outside  of  a  soil  specimen.   The  volume  of  any 
pores  that  are  completely  enclosed  by  surrounding  solids  cannot  be 
determined.   Moreover,  the  method  will  only  determine  the  volume  of 
intrudable  pores  that  have  an  apparent  diameter  corresponding  to  a 
pressure  within  the  pressuring  range  of  the  testing  instrument. 

Further  discussion  of  the  assumptions  and  sources  of  error  of  the 
mercury  intrusion  technique,  including  the  concern  of  whether  or  not 


11 


the  soil  structure  breaks  down  during  intrusion,  the  assumptions  of  con- 
stant contact  angle  and  surface  tension  values,  non-reaction  of  the 
mercury  with  the  soil,  and  incompressibility  of  the  soil  fabric,  among 
others,  is  found  in  Ritter  and  Drake  (19A5)  ,  Rootare  (1968),  Winslow 
(1969),  Diamond  (1970),  Bhasin  (1975),  Garcia-Bengochea,  et  al.  (1980), 
Kenney  (1980),  Reed,  et  al .  (1980),  and  a  proposed  ASTM  standard,  sub- 
conimittee  D18.14. 

2.1.3  The  Applications  of  Pore  Size  Distribution 
To  predict  the  permeability  of  porous  materials  is  one  of  the 
early  applications  of  the  PSD.   As  is  well  known,  seepage  occurs  only 
through  the  pores,  and  what  is  directly  related  with  the  permeability 
are  the  characteristics  of  pores  themselves.   Moreover,  permeability 
is  extremely  sensitive  to  subtle  changes  in  soil  structure  and  these 
changes  are  not  accurately  reflected  by  the  bulk  volumetric  pore  para- 
meters (e.g.,  the  porosity  or  void  ratio  of  the  sample).   This  led 
to  the  attempts  to  use  PSD  to  predict  the  permeability  of  porous 
material.   Among  those  are  Childs  and  Collis-George  (1950) ,  Marshall 
(1958),  Millington  and  Quirk  (1959),  Kunze,  et  al.  (1968),  Murota 
and  Sato  (1969),  Klock,  et  al.  (1969),  Green  and  Corey  (1971),  Jack- 
son (1972),  and  Garcia-Bengochea,  et  al.  (1979).   However,  only  the 
last  reference  clearly  demonstrated  the  usefulness  of  the  PSD  mea- 
surements on  interpreting  changes  in  soil  fabric  and  predicting  how 
these  changes  might  affect  the  permeability  of  soils. 

Another  application  of  PSD  is  to  predict  the  frost  heave  rate 
(or  amount)  in  soils.   References  include  Csathy  and  Townsend  (1962), 
Zoller  (1973),  Gaskin  and  Raymond  (1973),  and  Reed,  et  al.  (1979). 
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Among  other  applications,  Lohnes ,  Tuncer  and  Demirel  (1976)  re- 
lated PSD  to  weathering  and  mineralogy  of  loess  soils.   Ahmed,  et  al. 
(1974)  related  PSD  to  the  strength  of  clayey  soils. 

Finally,  based  on  PSD  studies  at  Purdue  University  during  the 
past  decade,  an  understanding  of  the  mechanism  of  compressibility 
and  that  of  liquefaction  of  sandy  soils  may  be  possible. 

2.2   PERMEAEILITY  MODEL 

2.1.1   Relation  of  Permeability  to  Pore  Geometry 

Since  permeability  is  a  characteristic  physical  property  of  a 
porous  medium,  it  would  seem  only  reasonable  to  assume  that  it  is 
related  in  some  functional  way  to  certain  measureable  properties  of 
the  soil  pore  geometry,  e.g.,  porosity,  pore  size  distribution,  internal 
surface  area,  etc.   However,  numerous  attempts  to  discover  a  functional 
relation  of  universal  applicability  have  so  far  met  with  disappointing 
results  (Hillel,  1980). 

Perhaps  the  simplest  approach  is  to  seek  a  correlation  between 
permeability  and  porosity  (or  void  ratio) .   However,  Scheidegger  (1957) 
and  Hillel  (1980),  among  others,  concluded  that  there  is  no  simple 
correlation  between  permeability  and  porosity  because  of  the  strong 
dependence  of  flow  rate  on  the  width,  continuity,  shape,  and  tortuo- 
sity of  the  conducting  channels.   Thus,  a  medium  composed  of  numerous 
small  pores  with  a  high  total  porosity  is  likely  to  exhibit  a  lower 
saturated  permeability  than  a  medium  composed  of  fewer  large  pores 
with  a  lesser  porosity. 

Prediction  of  permeability  from  grain  size  distribution  has  also 
been  used  for  a  long  time  (Hazen,  1911) .   Such  predictive  equations 
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have  had  only  limited  success  for  coarser  soils  and  have  proved  to  be 
unsatisfactory  for  fine-grained  soils.   Efforts  to  refine  this  approach 
by  introducing  empirical  grain  shape  and  packing  parameters  have  not 
won  general  acceptance  (Hillel,  1980)  . 

Numerous  theoretical  models  have  been  introduced  to  represent 
porous  media  by  a  set  of  relationships  that  are  amenable  to  mathemat- 
ical treatment.   Basically,  there  are  two  main  groups  of  those  models. 
The  first  is  based  on  a  generalization  of  Kozeny's  approach  (Carman, 
1S>36)  for  saturated  and  unsaturated  porous  media.   This  theory  is  based 
on  the  concept  of  a  hydraulic  radius,  i.e.,  a  characteristic  length 
parameter  presumed  to  be  linked  with  the  hypothetical  channels  to  which 
the  porous  medium  is  thought  to  be  equivalent.   The  derivation  and 
critique  of  this  theory  can  be  found  in  Scheidegger  (1957)  . 

The  second  group  employed  the  capillary  model  using  the  PSD  to 
calculate  permeability.   Scheidegger  (1957)  gave  a  comprehensive  re- 
view of  these  models,  including  the  straight  capillaric,  parallel, 
serial,  and  branching  models.   He  pointed  out  that,  in  general,  natural 
porous  media  are  extremely  disordered,  so  that  it  seems  a  rather  poor 
procedure  to  represent  them  by  something  which  is  intrinsically  ordered. 
He  therefore  suggested  that  the  preferred  model  of  a  porous  medium 
should  be  based  on  statistical  concepts. 

A  statistical  flow  model,  as  shown  in  Figure  2-1,  was  first  pre- 
sented by  Childs  and  Collis-George  (1950)  and  later  modified  by 
Marshall  (1958)  and  Millington  and  Quirk  (1959).   Conceptually,  their 
approach  rests  on  the  assumptions  that  (1)  a  soil  contains  distinct 
pores  of  various  sizes  which  are  randomly  distributed  in  space,  and 
(2)  when  adjacent  sections  of  the  soil  are  brought  into  contact,  the 
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Figure  2-1      Schematic  Representation  of  a 
Sectioned  Soil 
(After  Hillei,  1980) 
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overall  quantity  of  flow  across  the  section  depends  statistically  on 

the  number  of  pairs  of  interconnected  pores  and  geometrically  on  their 

configurations.   The  quantity  of  flow  of  each  pair  of  interconnected 

pores  is  determined  by  the  narrower  of  the  two  pores. 

The  original  equation  by  Childs  and  Collis-George  (1950)  was 

s=R  6=R 
k  =  c  -2-^   Z   S  5  f (s)  ds  f(6)  d6  (2-9) 

s=o  6=o 

where        k  =  permeability 

c  =  a  matching  constant 
p  =  density  of  water 
g  =  gravitational  acceleration 

M  =  coefficient  of  absolute  viscosity  of  water 
f(s)  ds  =  the  partial  pore  space  occupied  by  pores  with  radii 

of  s  ->  s  +  ds 
f(6)  dS    =   the  partial  pore  space  of  pores  with  radii  of 
6  -»■  6  +  d6 

R  =  maximum  pore  radius 
Note  that  although  theoretically  the  calculation  based  on 
Poiseuille's  law  does  not  require  a  matching  factor,  in  actual  prac- 
tice a  matching  factor  is  needed  to  adjust  the  computed  and  measured 
permeability  at  full  saturation.   The  method  of  calculation  was  sub- 
sequently improved  by  Marshall  (1958),  Millington  and  Quirk  (1959), 
Kunze  et  al.  (1968),  Green  and  Corey  (1971),  and  Jackson  (1972). 
Basically  these  developments  are  different  only  by  some  minor  arbi- 
trary factors  and  the  ways  of  calculation.   Moreover,  Garcia-Bengochea 
(1978)  presented  a  somewhat  different  equation  based  on  regression 
analysis  of  his  data  as  follows: 
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k  =  c   •  PSP^  (2-lOa) 

s 

where       b  =  parameter  of  regression 

c  =  a  constant,  called  shape  factor 
s 

n  n  „ 
and  PSP=  Z  I  d  f(d.)  f(d.)  (2-lOb) 

i  j       ^     -^ 

d  =  smaller  of  d.  and  d. 
fCd,),  f(d.)  =  volumetric  frequency  of  occurrence  of  pores  with 
diameter  of  d  ->  d  +  Ad  ;  conceptually  it  equals 
f(s)ds  or  f(6)d6  in  Equation  (2-9). 
As  noted  by  Garcia-Bengochea,  et  al.  (1979),  as  b  ->-  1,  Equation 
(2-10)  becomes  the  Child  and  Collis-George's  (1950)  model.   Since  ex- 
perimentally b  ranges  from  1.67  to  A. 95,  Equation  (2-10)  is  indeed  an 
empirical  equation. 

2.2.2  Proposed  Permeability  Model 

From  the  previous  discussion,  it  is  clear  that  Child  and  Collis- 
George's  (1950)  model  and  its  modifications  are  more  generally  appli- 
cable than  those  based  on  earlier  models.   However,  an  arbitrary  match- 
ing factor  or  regression  parameter,  as  appears  in  Equation  (2-9)  and 
Equation  (2-10)  respectively,  is  necessary  to  reduce  or  to  eliminate 
the  deviations  between  predicted  and  measured  permeability.   One  of 
the  objectives  of  the  present  study  is  to  propose  a  modified  model 
which  eliminates  the  deviation  without  the  necessity  of  an  arbitrary 
factor  or  parameter. 

Considering  a  homogeneous  soil  having  interconnected  pores  which 
are  randomly  distributed  in  space,  the  probability  of  pores  of  diameter 
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X  -»■  X  +  dx  is  f(x)  dx,  where  f(x)  is  the  pore  size  density  function  of 
the  soil.   Note  that  no  assumption  is  needed  for  f (x) .   However,  in  the 
present  study,  f(x)  is  defined  by  Equation  (2-7a).   The  term  f(x)  dx, 
represents  the  ratio  between  the  pore  volume  occupied  by  pores  of  dia- 
meter X  ->-  X  +  dx  and  the  total  pore  volume.   Note  that  if  a  cross  sec- 
tion through  soil  mass  is  considered,  the  above  volumetric  ratio  is 
equal  to  the  areal  ratio  between  the  pore  area  of  diameter  x  ^  x  +  dx 
and  the  total  pore  area. 

Now,  consider  a  soil  column  with  a  thickness  Ay  and  two  cross 
sections  with  an  identical  pore  size  density  function  (see  Figure  2-2). 

The  probability  of  pores  of  diameter  x.  -*-  x.  +  dx .  on  cross  section  i 
^         ^     ^  111 

being  connected  to  pores  of  diameter  x.  ^  x.  +  dx.  on  cross  section  j, 
P(x.,  X.),  has  the  following  extreme  cases: 

(1)  For  Ay  >>  x,  the  connections  between  pores  on  the  two  cross 
sections  can  be  reasonably  assumed  to  be  completely  random. 
Hence 

P(x.,  X.)  =  f(x^)  f(x.)  dx^  dx.  (2-11) 

(2)  For  Ay  -»-  0,  the  connections  between  pores  considered  are  com- 
pletely correlated. 

Hence 

P(x.,  X.)  =  f(x.)  dx.  or  f(x.)  dx.  (2-12) 

1   J       i'    1       J    J 

For  case  1,  following  the  similar  derivation  by  Childs  and 
Collis-George  (1950)  or  Garcia-Bengochea  (1978)  ,  we  have 

k  =  -^  /    /     ^^  f(xJ  f(x.)  dx.  dx.  (2-13) 

^^^   Jo        Jo  i     J     1    J 
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where  y   =  unit  weight  of  water 

\i  -   coefficient  of  absolute  viscosity  of  water 

n  =  porosity  of  soil 

X  =  equivalent  diameter  of  capillary  element. 

Equation  (2-13)  is  essentially  equal  to  Equation  (2-9)  except  that  no 

matching  factor  is  provided.   This  equation  is  also  called  the  Marshall 

model  by  Garcia-Bengochea  (1978). 

For  case  2 ,  in  a  similar  manner,  we  have 

2  /"^ 
k  =  ^  /    x^  f(x)  dx  (2-lA) 

Jo 

Equation  (2-14)  is  called  the  capillary  model  by  Garcia-Bengochea 
(1978). 

Since  we  are  concerned  with  the  effect  of  pore  geometry  changes 
on  the  permeability,  it  is  more  relevant  to  consider  Ay  to  be  of  the 
same  order  of  magnitude  as  the  pore  diameter.  Therefore  the  connec- 
tions between  pores  on  the  two  cross-sections  are  partially  correlated. 

A  reasonable  assumption  to  account  for  the  partially  correlated 

connections  is  that  pores  of  diameter  v;  -»-  x  +  dx   on  cross-section  i 

i    i     i 

have  a  greater  probability  of  being  connected  to  pores  of  the  same  dia- 
meter (i.e.,  ^i  "*"  ^i  ^   i)  on  cross-section  j,  than  that  of  being 
connected  to  the  other  different-sized  pores  on  cross-section  j . 
Based  on  this  idea,  a  connection  rule  is  illustrated  in  Figure  2-3. 

The  probability  of  pores  of  diameter  x.  "*"  ^^  +  dx^  on  cross- 
section  i  being  connected  to  pores  of  diameter  x.  ->-  x .  +  dx.  on  cross- 
section  j  is 

P(^.  X  )  =  g(y,  X    X  ).  •  f(x.)  •  f(x.)  dx,  dx.  (2-15) 

J  ■'-J       1       jij 
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and 


2  /*°°  /"^^ 
^/   /   x2  .  g(y  ,  X..  xj)-f(x^)  '  f(x.)  dx^  dx.       (2-16) 


/O   JO 


where  g(y,  x  ,  x.)  =  connecting  function  and 

y  =  the  third  directional  variable  (length) . 
If  we  also  consider  a  correction  factor,  t(y,  x.,  x.)  <  1,  to 
account  for  the  tortuosity  of  flow,  then  the  permeability  becomes 

k  =  ^/    /    x-g-t-  f(x^)  .  f(xj  dx.  dx_.  (2-17) 


■'o     Jo 


Based  on  the  above  development,  we  now  introduce  a  governing 
function  G(x.)  to  account  for  g  and  t  as  follows  (see  Figure  2-3) 


;(x.)=  e  -  (  "'  V  } 


(2-18) 


where  V  =  V(x.) 

=  an  undetermined  function  of  pore  size  to  account  for  the 
third  directional  variable  and  the  tortuosity  of  flow. 
Note  that  pore  size  density  functions  f(x.)  and  f(x.)  are  identical. 
And  random  variables  X.  and  X.  are  symmetrical  in  Equation  (2-18) . 

For  any  pores  of  diameter  x.  -»■  x.  +  dx,,  the  governing  function 
is  bell-shaped  and  symmetrical  with  a  mean  value  of  x.  -^  x .  +  dx., 
as  shown  in  Figure  2-3.   Therefore  Equation  (2-18)  is  a  family  func- 
tion and  each  member  of  the  family  is  determined  by  a  pair  of  para- 
meters, i.e.,  X  and  V. 

Based  on  the  above  discussion,  Equation  (2-17)  can  be  refined  as 
follows : 

2  /oo  ^  CO 

^'^\     ]         5  •  G^^j)  •  f(Xi)  •  f(Xj)  dx^  dxj  (2-19) 
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Equation  (2-19)  is  the  proposed  permeability  model  which  has  the 
following  characteristics: 

(1)  a  probabilistic  basis; 

(2)  no  arbitrary  matching  factor  required;  and 

(3)  a  tortuosity  factor  being  considered. 
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CHAPTER  III   TESTING  PROGRAM  AND  EXPERIMENTAL  PROCEDURES 

3.1   TESTING  PROGRAM 

The  testing  program  was  designed  to  accomplish  the  objectives 
mentioned  in  Section  1.2.   The  flow  chart  of  the  program  is  shown  in 
Figure  3-1.   The  soils  were  selected  in  order  to  extend  the  results  of 
the  previous  permeability-pore  size  distribution  study  at  Purdue  Univer- 
sity by  Garcia-Bengochea  (1978),  i.e., to  extend  the  prediction  model 
to  wider  ranges  of  permeability  and  to  characterize  different  fabrics 
of  sandy  soils.   The  compaction  variables  considered  in  this  study  are 
water  content,  energy  level,  and  compaction  method. 

The  pore  size  distributions  (PSD)  of  the  selected  soils  were  de- 
termined by  the  mercury  intrusion  method.   For  clay-sand  mixtures, 
following  the  permeability  test,  the  PSD  specimens  were  trimmed  and 
dried  by  a  f reeze-drying  technique  and  were  made  ready  for  the  PSD 
test.   For  pure  sand,  the  problem  of  determining  its  pore  size  distri- 
bution is  in  sampling,  i.e., how  to  get  a  PSD  specimen  without  destroy- 
ing its  fabric.   In  the  present  study,  this  problem  has  been  solved  by 
using  a  resin  as  a  bonding  material.   A  small  amount,  0.5%  by  weight, 
of  Varcum  1364  phenolic  resin  made  by  Reichhold  Chemicals,  Inc.  was  mixed 
with  sand  before  compaction.   After  compaction  the  entire  sample  and 
mold  were  placed  in  the  oven  at  a  temperature  of  170  C  for  1  h. 
The  sample  was  then  removed  and  placed  in  a  desiccator  several  hours 
until  it  cooled.   The  oven-dried,  compacted  samples  were  found  to  have 
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Figure  3-1    Flow  Chart  of  the  Testing  Program 
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about  5%  shrinkage  in  volume  but  no  change  in  total  weight.  Wins low( 1980) 
stated  that  no  chemical  reaction  between  sand  and  powdered  phenolic  resin 
occurs  during  oven-drying  because  the  resin  only  acts  as  a  glue.  Moreover, 
the  results  of  the  PSD  tests  on  specimens  trimmed  from  these  oven-dried 
samples  are  reproducible.   Therefore  this  technique  is  believed  to  be 
the  best  way  to  obtain  PSD  specimens  of  pure  sand.   During  some  pre- 
liminary work,  the  minimum  amount  of  resin  which  provides  sufficient 
bonding  and  workability  was  found  to  be  0.5%  by  weight  of  sand. 

The  following  sections  desi^ribe  the  soils  studied,  compaction 
methods,  permeability  tests,  specimen  drying  methods,  and  the  PSD  tests. 

3.2  SOILS  STUDIED 

Sand  and  sand-clay  mixtures  used  in  this  study  have  been  arbi- 
trarily called  Soils  A,  B,  and  C,  respectively.   The  specifications 
of  these  three  soils  are  given  in  Table  3-1.   The  clay  used  is  Edgar 
plastic  kaolin  with  the  following  index  properties  (Garcia-Bengochea, 
1978):   specific  gravity  of  solids  =  2.65,  liquid  limit  =  59.',  plas- 
ticity index  =  22%,  and  the  clay  fraction  (%  less  than  2  ym)  =  76%. 
The  sand  used  is  a  uniform  Ottawa  sand.   The  grain  size  distribution 
curves  and  the  classification  properties  of  Soils  A,  B,  and  C  are 
shown  in  Figure  3-2  and  Table  3-1  respectively. 

3.3  COMPACTION 

3.3.1   Pluvial  Compaction  for  Soil  A 

For  Soil  A,  0.5%  powdered  phenolic  resin  was  mixed  with  Ottawa 
sand   before  compaction.   Pluviation  was  used  to  compact  the  sample 
in  a  standard  Proctor  compaction  mold  (ASTM  D698)  with  a  diameter  of 
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102  mm  and  116  mm  high.   Figure  3-3  illustrates  the  scheme  of  the  plu- 
viation.   The  sand-resin  mixture  was  poured  from  a  1000  ml  flask 
through  a  nozzle  made  from  rubber  stopper  containing  one  hole  of  either 
6.8  or  3.8  mm  diameter.   The  sand-resin  mixture  was  free  to  fall  to 
the  surface  of  the  sample  while  the  flask  was  rotated  at  a  rate  of 
approximately  15  rpm.   The  density  of  the  sample  obtained  was  a  func- 
tion of  the  intensity  of  particle  rain,  which,  in  turn,  depended  on 
hole  size,  and  the  dropping  height.  In  the  present  study,  the  sample  was 
made  in  five  layers  with  height  of  drop  for  each  layer  being  108  mm. 
In  order  to  keep  the  same  height  of  drop,  the  position  of  flask  was 
adjusted  for  each  layer.   The  choice  of  dropping  height  and  hole  size 
depended  on  the  density  of  the  compacted  sample  desired.   For  this 
study,  they  were  selected  by  trial  and  error. 

3.3.2  Kneading  Compaction  for  Soils  B  and  C 

An  electrically  driven  semi-automatic  kneading  compactor  shown 
in  Figure  3-4  was  used  to  compact  the  samples  in  a  standard  Proctor 
compaction  mold.   The  compaction  mold  was  also  used  as  the  permeameter 
cell  in  order  to  reduce  the  possibility  of  seepage  along  the  walls  of 
the  mold.   The  kneading  compaction  procedures  described  by  Garcia- 
Bengochea  (1978)  were  used  in  this  study  with  the  exception  of  the 
definition  of  the  energy  level.  In  the  present  study,  the  density  level 
was  defined  by  the  standard  and  modified  Proctor  compaction  as  speci- 
fied in  ASTM  D698  and  D1557.  The  gage  pressure  set  on  the  keneading 
compactor  at  each  water  content  was  selected  by  trial  and  error  in  an 
attempt  to  obtain  the  same  density  as  indicated  on  the  Proctor  compac- 
tion curves.   In  Garcia-Bengochea's  (1978)  study,  the  gage  pressure 
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Figure  3-3       Schematic  Diagram  of  Pluviation  System 


30 


Figure  3-4   Kneading  Compactor 
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was  set  to  the  same  value  for  each  water  content  to  define  an  energy 
level,  as  shown  in  Figure  3-5.   The  reasons  for  the  present  approach 
are  (1)  the  density-water  content  curves  for  different  compaction 
methods  may  not  have  the  same  shape,  and  (2)  the  Proctor  compaction 
curve  is  more  common  in  soils  engineering  practice  for  specifying  com- 
paction energy.   Following  the  present  procedure,  calibration  curves  for 
converting  gage  pressure  to  the  actual  foot  pressure  are  not  necessary. 
Hence  the  problem  of  variation  of  calibration  curves  by  seasonal  changes 
of  temperature  and  humidity,  as  noted  by  Garcia-Bengochea  (1978),  has 
been  eliminated. 

After  mixing  and  a  one  day  curing  period,  the  soil  was  placed  in 
the  Proctor  mold  and  compacted  in  five  layers  of  approximately  equal 
thickness.   Thirty  tamps  of  the  compaction  foot  were  applied  for  each 
layer.   Prior  to  compaction  of  each  layer  except  the  first  one,  the 
previous  layer  was  scarified,  and  then  fresh  soil  was  added.   The 
scarification  of  the  layers  prevented  separation  between  the  layers. 
Once  the  compaction  procedure  was  completed,  the  collar  of  the  mold 
was  removed  and  the  soil  sample  was  trimmed  flush  with  the  top  of  the 
Proctor  mold. 

3.4  PERMEABILITY  TESTS 

Basically,  the  permeability  testing  equipment  and  procedures  used 
in  this  study  followed  Garcia-Bengochea  (1978). 

For  soil  A,  the  compacted  sample  in  the  Proctor  mold  was  connected 
to  the  saturation  and  permeation  device  after  oven  drying.   For  soils 
B  and  C,  the  sample  in  the  mold  was  connected  to  the  saturation  and 
permeation  device  directly  after  completion  of  compaction  process. 
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Figure    3-5       Compaction  Curves  for 90%  Silt  and  10% 
Kaolin 
(After  Garcia-Bengocheo ,  1978) 
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The  saturation  and  permeation  device  is  a  falling  head  permeameter 
system  with  provisions  for  applying  both  vacuum  and  back  pressure.   The 
apparatus  shown  in  Figure  3-6  consists  of  three  air  pressure  regulators, 
one  air  pressure  gage,  two  standpipes  and  air  pressure  and  vacuum  lines. 
The  two  permeameters  are  nearly  identical  except  for  the  inside  diameters 
of  the  standpipes.   The  permeameter  used  for  measuring  the  more  permeable 
soils  has  12.7  mm  inside  diameter  (I.D.)  standpipes,  while  the  other  has 
6,3  mm  I.D.  standpipes.   The  standpipes  can  be  pressurized  simultane- 
ously from  one  or  independently  from  two  regulators.   Meter  sticks  were 
used  to  record  the  water  elevation  heads.   The  device  is  designed  for 
both  saturation  and  permeation.  A  permeameter  cell  is  shown  in  Figure  3-6. 

For  soil  C,  a  vacuum  of  about  34.5  kPa  (5  psi)  absolute  pressure 
was  placed  on  the  top  of  the  sample  while  the  bottom  end  was  under  approx- 
imately 1  m  of  water  head.   This  condition  was  maintained  for  about 
10  hours  which  allowed  de-aired  water  to  flow  upward  through  the  sample 
and  allowed  air  bubbles  to  escape.   For  soils  A  and  B,  due  to  their 
relatively  high  permeabilities,  the  above  vacuum-saturation  process 
was  accomplished  within  an  hour. 

After  vacuum  saturation,  back  pressure  was  applied  in  increments 
of  34.5  kPa  (5  psi)  to  both  ends  of  the  sample.   Each  increment  was 
maintained  for  10  minutes,   ^^fhen  the  back  pressure  reached  345  kPa 
(50  psi),  it  was  maintained  for  about  15  hours. 

In  the  above  saturation  process,  no  B- check  is  possible  to  assure 
that  the  sample  is  fully  saturated.   However,  as  noted  by  Garcia- 
Bengochea  (1978)  ,  the  vacuum  saturation  consistently  resulted  in  a  sat- 
uration level  of  at  least  95%  for  the  clay-silt  mixtures.   There  is  no 
measurement  of  the  degree  of  saturation  after  vacuum  saturation.   Accord- 


34 


Figure  3-6   Saturation  and  Permeation  Device 
and  Permeameter  Cells 
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Ing  to  Black  and  Lee  (1973)  ,  the  back  pressure  required  to  reach  a  99':: 
saturation  is  about  300  kPa  (44  psi)  for  initial  degree  of  saturation 
of  95%.   The  corresponding  time  required  to  reach  equilibrium  is  about 
4  hours.   For  the  saturation  procedure  described  above,  it  is  believed 
the  degree  of  saturation  is  at  least  99%,  which  is  appropriate  for 
the  practical  purposes  of  permeability  measurement.   To  check  this 
inference,  a  sample  was  tested  following  the  above  procedures  except 
that  a  back  pressure  of  483  kPa  (70  psi)  was  applied  during  the  satur- 
ation and  permeation  phases.   No  measureable  difference  was  found 
between  the  measured  permeabilities. 

At  the  completion  of  the  saturation  procedure,  the  falling  head 
permeability  tests  were  conducted.   The  direction  of  flow  was  from 
top  to  bottom  across  the  sample.   A  initial  head  difference  of  about 
100  cm  or  45  cm  was  created  between  the  two  standpipes  while  a  345  kPa 
(50  psi)  back  pressure  was  maintained  on  both  standpipes.   Readings 
of  the  elevation  heads  in  the  water  columns,  time,  and  air  temperature 
was  recorded. 

Following  permeation,  the  sample  was  dissected  to  check  for  uni- 
fomity  and  trimmed  for  freeze  drying  and  pore  size  distribution  tests. 
As  noted  by  Garcia-Bengochea  (1978),  the  back  pressure  on  the  sample 
should  not  be  released  too  rapidly  in  order  to  prevent  possible  sample 
disturbance.   Thus  each  decrement  of  34.5  kPa  (5  psi)  was  maintained 
for  15  minutes.   When  atmospheric  pressure  was  reached,  the  sample  was 
left  in  the  mold  for  at  least  2  hours  before  dissection  and  specimen 
trimming. 
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3.5   SPECIMEN  PREPARATION  AND  DRYING  AND  PSD  TESTS 

For  Soil  A,  specimens  for  the  PSD  test  were  trimmed  from  oven- 
dried  compacted  samples.   The  dimensions  of  Individual  PSD  specimens 

were  approximately  12  x  8  x  8  mm.   After  trimming,  the  PSD  specimens 

were  placed  in  an  oven  for  ten  minutes  in  order  to  get  rid  of  moisture 

which  may  have  been  absorbed  by  the  specimens  during  trimming.   The  PSD 

specimens  were  then  kept  in  glass  bottles  and  placed  in  the  desiccator. 

For  Soils  B  and  C,  after  permeability  testing,  the  PSD  specimens 
were  carefully  trimmed  with  a  razor  blade  into  cubes  approxin.ately 
8  mm  on  a  side.   The  method  for  drying  Soils  B  and  C  in  this  study  was 
the  freeze  drying  technique.   The  method  consists  of  quickly  freezing 
a  small  soil  specimen  to  cryogenic  temperatures,  followed  by  vacuum 
drying  the  frozen  sample  to  remove  water  by  sublimation.   This  process 
eliminates  the  air-water  menisci  shrinkage  effect,  which  may  alter  the 
soil  structure,  by  replacing  it  with  an  air-solid  system.   The  small 
shrinkage  that  does  occur  (less  than  5%  as  determined  by  testing  a  few 
specimens  initially)  after  f reeze-drying  indicates  that  minimal  sample 
disturbance  probably  occurred. 

The  freeze-drying  apparatus  and  procedures  used  in  this  study 
followed  Garcia-Bengochea  (1978) .   Appendix  A-1  provides  details  of 
this  apparatus  and  procedures. 

After  the  specimens  were  dried,  the  PSD  tests  were  conducted.   The 
apparatus  and  procedures  used  in  this  study  was  similar  to  those  used  by 
Garcia-Bengochea  (19  78)  and  White  (1980).   Appendix  A-2  provides  de- 
tails of  PSD  testing  equipment  and  step  by  step  procedures.   However, 
it  is  worthwhile  to  point  out  here  that  both  penetrometers  (with 
intrusion  capacity  of  0.2  cc  and  1.2  cc  respectively)  were  used  in 
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the  present  study.   This  is  because  the  small  penetrometer  (with  capa- 
city 0.2  cc) ,  as  used  by  Garcia-Bengochea  (1978)  and  White  (1980), 

cannot  properly  measure  the  PSD  of  the  high  permeability  soils  A  and 
B.   Also  note  that  the  large  penetrometer  (with  capacity  1.2  cc)  can 
not  be  placed  in  the  Amico  porosimeter;  therefore  it  can  only  go  through 
low  pressure  intrusion  which  measures  the  distribution  of  the  large 
size  pores  (20-  500  ym) .   For  these  reasons,  a  combination  of  data 
measured  by  both  penetrometers  is  used  to  show  the  distribution  of 
the  entire  range  of  pore  sizes. 
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CHAPTER  IV  PRESENTATION  AND  DISCUSSION 

4.1   COMPACTION 

4.1.1  Sample  Nomenclature 

Each  sample  tested  in  this  study  was  assigned  a  code  name  based 
on  the  soil  type,  compactive  effort,  compaction  method,  and  relative 
water  content  with  respect  to  its  optimum.  This  code  designation  is 
listed  in  Table  4-1. 

4.1.2  Dry  Density  versus  Water  Content  Relationships  for  Soils 
B  and  C. 

The  dry  density  versus  water  content  relationships  utilized  in 
this  study  were  determined  from  the  standard  and  modified  Proctor  com- 
paction tests.  The  specifications  for  these  impact  compaction  tests 
are  listed  in  Table  4-2.  The  Proctor  test  results  are  listed  in  Table 
4-3.  The  statistical  relationships  between  dry  density  and  water  con- 
tent obtained  from  a  least  squares  analysis  of  these  data  are  shown  in 
Figure  4-1.  Satisfactory  fit  to  the  data  is  indicated  by  the  high  co- 
efficients of  multiple  determination  listed  in  Table  4-3. 

A  California  type  kneading  compactor  was  used  to  compact  most  of 
the  samples  tested  (i.e..  Soils  B  and  C)  in  this  study.   The  kneading 
compaction  dry  density  and  water  content  points  used  in  the  testing 
program  are  shown  in  Figure  4-1.   The  gauge  pressures  utilized  for 
each  of  these  points  are  listed  in  Table  4-4.   For  Soil  C,  both  re- 
sults of  high  and  medium  compactive  efforts  are  presented.   However, 
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TABLE  4-1   Sample  Code  Designation 

Soil  Type 

0  Soil  A  (Ottawa  sand  only) 

SCI  Soil  B  (90%  Ottawa  sand  +  10%  Kaolin) 

SC3  Soil  C  (70%  Ottawa  sand  +  30%  Kaolin) 

Compact ive  Effort 

For  Soil  A: 

L  Relatively  low  effort 

H  Relatively  high  effort 

For  Soils  B  and  C: 

M  Medium  effort  (specified  in  ASTM  D698) 

H  High  effort  (specified  in  ASTM  D1557) 

Water  Content 

0  Optimum  water  content 

0-D  Very  slightly  dry  of  optimum 

D  Dry  of  optimum 

W  Wet  of  optimum 

Compaction  Method 

nil  Kneading  compaction 

P  Pluvial  compaction 

Example: 

SCIMO:   Soil  B  (90%  Ottawa  sand  +  10%  Kaolin) ,  compacted  by 
medium  effort  at  optimum  water  content  by  kneading 
compaction. 
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LEGEND 
I.  Statistical  curve,  std.  Proctor  comp.' 
2.  Statistical  curve,  mod.  Proctor  comp." 
3.2ero-Qir-voids  curve  (100%      saturntron) 
A   Sample  point  by  kneading   comp.,  Soil  B 
G    Sample  point  by  kneading  comp..  Soil  C 
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Figure  4-1         Dry  Density  versus  Water  Content  Curves 
for  Soils  B  and  C 


Curves  are  somewhat  displaced  and  different  in  shape  from  what 
would  have  been  expected,  but  are  based  upon  real  data. 


TABLE  4-4  Kneading  Compaction  Variables  for  Each  Sample  Point  of 
Soils  B  and  C 
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Sample   Sample 
Code   Point 


Water  Content   Gauge  Pressure 
(%)  kPa  (psi) 


Achieved 
Dry  Density 

(Mg/m^) 


1 

9.1 

55.2  (8) 

1.87 

SC3MD 

2 

9.2 

55.2 

1.87 

3 

8.4 

55.2 

1.81 

1 

11.0 

.  55.2 

1.96 

SC3M0 

2 

10.9 

55.2 

1.96 

1 

12.5 

55.2 

1.90 

SC3MW 

2 

12.5 

55.2 

1.91 

1 

6.5 

193.2  (28) 

1.95 

SC3HD 

2 

6.6 

193.2 

1.95 

3 

6.2 

193.2 

1.94 

1 

8.0 

207.0  (30) 

1.99 

SC3H0-D 

2 

7.9 

207.0 

1.99 

1_ 

6.2 

34.5  (5) 

1.88 

SCIMD 

2 

6.1 

34.5 

1.88 

3 

5.9 

34.5 

1.87 

1 

8.0 

41.4  (6) 

1.89 

SCIMO 

2 

8.1 

41.4 

1.89 

3 

8.0 

41.4 

1.87 

1 

11.0 

34.5 

1.89 

SCIMW 

2 

11.1 

34.5 

1.88 
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for  Soil  B,  only  the  result  of  medium  compactive  effect  of  this  soil 
was  presented  because  1)  the  compaction  result  was  little  affected  by 
the  magnitude  of  the  impact  compaction  effort;  and  2)  the  appropriate 
gauge  pressures  necessary  to  fit  the  results  of  kneading  compaction 
to  the  desired  density  as  defined  by  the  modified  Proctor  compaction 
curve  were  not  able  to  be  found  in  the  present  study. 

4.1.3   Pluvial  Compaction  Results  for  Soil  A 

The  results  of  Soil  A  obtained  from  the  pluvial  compaction  are 
listed  in  Table  4-5.   The  results  indicate  that  the  desired  density 
can  be  replicated  by  controlling  the  diameter  of  hole  and  the  height 
of  drop. 

4.2   PERMEABILITY 

The  permeability  was  measured  by  a  falling  head  apparatus  under  a 
back  pressure  of  345  kPa  (50  psi) .   The  hydraulic  gradient  during  the 
tests  ranged  between  2  to  5.   The  permeabilities  were  calculated  from 
the  usual  falling  head  equation  (Taylor,  1948) 

t 

where  k  =  permeability 

a  =  cross-sectional  area  of  the  standpipes 

A  =  cross-sectional  area  of  sample 

At  =  time  increment 

h  =  total  head  at  zero  time 
o 

h  =  total  head  after  time  At. 
Figure  4-2  shows  a  schematic  diagram  of  the  falling  head  test. 
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TABLE  4-5  Pluvial  Compaction  Variables  for  Each  Sample  Point  of  Soil  A 


Compacted 

Sample 

Sample 

Diameter  of 

Height  of 

Code 

Point 

Hole  (mm) 

Drop  (mm) 

Number     Achieved 
of  Layers   Dry  Density 

(Mg/m-^) 


OPL 


1 
2 

3 
4 
5 
6 


6.8 
6.8 
6.8 
6.8 
6.8 
6.8 


108 

5 

1.62 

108 

5 

1.62 

108 

5 

1.62 

108 

5 

1.62 

108 

5 

1.62 

108 

5 

1.62 

OPH 


1 
2 
3 
4 
5 
6 


3.8 
3.8 
3.8 
3.8 
3.8 
3.8 


108 

5 

1.72 

108 

5 

1.72 

108 

5 

1.72 

108 

5 

1.72 

108 

5 

1.72 

108 

5 

1.72 
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Figure  4-2       Schematic  Diagram  of  Falling  Head  Permeability 
Test  (After  Garcia-Bengochea,  1978) 
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It  is  evident  from  the  figure  that 


h   =  (h,   +  u  )  -  (h„   +  u  ) 
o      Ao     o       Bo     o 

=  h^   -  h„  .  (4-2) 

Ao     Bo 

where  u   =  back  pressure  on  the  svstem 
o 

and 


\   =  V  -  ^Bt  ^'-'' 

The  permeability  is  then  calculated  from  Equation  (4-1)  . 

The  variation  of  the  room  temperature  during  permeation  of  most 
of  the  permeability  tests  was  negligible  (e.g.,  +  0.4  C) ,  although  in 
a  few  tests  the  variation  was  up  to  +  1  C  or  so.   For  each  test,  an 
average  temperature  was  taken  and  the  following  equation  suggested  by 
Taylor  (1948)  was  used  to  calculate  and  to  report  the  permeability: 

^20°C  =  \  •  -^  (^-^) 

>^20°C 

where  ^^r^o^   =  permeability  at  temperature  20  C. 

k,  =  permeability  at  temperature  T  of  test 

y„  =  viscosity  of  water  at  temperature  T 

y^„o   =  viscosity  of  water  at  temperature  20  C. 

A  summary  of  the  permeability  test  results  of  the  compacted  sam- 
ples listed  in  Table  4-4  and  Table  4-5  is  given  in  Table  4-b .   In  this 
table  the  permeability  is  normalized  to  an  ASTM  standard  temperature 
of  20°C 

For  Soil  C,  the  dry-side  samples  had  permeab Llitics  one  to  two 
orders  of  magnitude  greater  than  the  optimum  and  wet-side  samples. 
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TABLE  4-6   Permeabilities  of  Samples  Tested 


Soil     Sample        Point  ^■yo°c  Average  k  o 

Type      Code         No.  (cm/sec)  (cm/sec) 

1  4.6  X  10~^ 

SC3MD  2  4.6  X  10"^  4.5  x  10~^ 


SC3M0  '       ,  .^  -  ^v.  7.3  x  10  ^ 


A 


SC3MW 


3 

4.2   X   10   ^ 

1 

7.1  X   lO"^ 

2 

7.5  X   O"'' 

1 

6.1  X   lo"'' 

2 

5.8   X    10~^ 

1 

1.8  X   10~^ 

3 

3.9 

X 

10-^ 

1 

1.3 

X 

10-^ 

2 

1.3 

X 

10-^ 

3 

1.2 

X 

10-^ 

1 

5.6 

X 

10-^ 

2 

5.0 

X 

10-^ 

5 

2.7 

X 

10-2 

6 

2.3 

X 

10-2 

5 

1.4 

X 

10-2 

6 

1.3 

X 

10-2 

6.0  X  lO"'' 


SC3HD  2  1.7x10^  1.9x10^ 

3  2.1  X  IQ-^ 

nrnun    n         ^  6 . 1  X  10  .  ,     ,^-7 

SC3H0-D  _  6 . 1  X  10 

2  6.0  X  10 

1  3.6  X  10-"^ 

SCIMD  2  3.0  X  lO"^  3.5  x  lo"^ 


SCIMO  2       1.3  X  10  "^  1.3  x  lo""^ 


SCIMW  *       "'"  "  ""  ,  5.3  X  10  '^ 


GPL  -'       ^.,  ^  ^v.  2.5  X  10  2 


OPH  '        ^.^  ^  .^  1,4  X  10"2 
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There  was  little  difference  between  the  permeabilities  of  the  optimum 
and  wet-side  samples.   These  trends  were  found  by  the  previous  research- 
ers, e.g.,  Lambe  (1954),  Mitchell,  Hooper  and  Campanella  (1965),  and 
Garcia-Bengochea  (1978)  .   However,  for  Soil  B  which  contains  only  10% 
clay  by  weight,  the  above  trend  is  not  very  clear;  i.e.,  the  difference 
between  the  dry-side  and  optimum  samples  was  small.   Two  possible 
reasons  are  1)  water  content  is  not  a  sensitive  variable  in  the  compac- 
tion of  this  soil,  as  indicated  in  Figure  4-1,  and  2)  the  dry-side 
samples  tested  are  not  suf  f  icie.\tly  dry  of  optimum  (see  Figure  4-1)  . 
Strictly  speaking,  there  is  no  clear  optimum  water  content  on  the 
standard  Proctor  compaction  curve  of  this  soil.   This  may  be  the  nature 
of  this  poorly  graded  sandy  soil. 

4.3   PORE  SIZE  DISTRIBUTION 

4.3.1  Data  Reduction 

As  mentioned  in  Section  2.1.1,  the  absolute  applied  pressures  can 
be  converted  to  apparent  intruded  pore  diameters  based  on  Equation  (2-1) 
This  step  requires  that  the  surface  tension  and  contact  angle  be  known. 
The  surface  tension  used  in  this  study  is  assumed  to  be  484  dynes/cm 
for  all  soils  used  (Kemball,  1946  and  Bhasin,  1975) .   The  contact  angle 
of  mercury  on  Kaolin  is  147   (Diamond,  19  70)  and  that  of  mercury  on 
quartz  (sand)  is  also  147   (Ellison,  Klemm  and  Schwartz,  1967).   Based 
on  these  reported  values  (both  at  25  C) ,  the  contact  angle  of  mercury 
on  Soils  B  and  C  is  assumed  to  be  147  .   For  Soil  A,  the  contact  angle 
of  mercury  on  phenolic  resin  was  determined  to  be  154   following  the 
method  suggested  by  Winslow  (1969,  1980).   Since  the  surface  tension 
of  phenolic  resin  is  much  lower  than  that  of  sand,  therefore  the  con- 
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tact  angle  of  mercury  on  Soil  A  was  assumed  to  be  154  . 

For  each  PSD  test,  the  applied  pressures  and  intrusion  data  were 
recorded.   The  pressures  recorded  during  the  filling  operation  were 
absolute  pressures,  while  those  recorded  during  the  porosimcter  oper- 
ation (called  high  pressure  intrusion)  were  gage  pressures.   As  noted 
by  Bhasin  (1975),  a  net  correction  of  76  kPa  (11  psi)  accounting  for 
atmospheric  pressure  and  negative  head  of  mercury  is  added  to  the  gauge 
pressures  to  obtain  the  absolute  pressures.   The  absolute  pressures 
are  next  converted  to  apparent  intruded  pore  diameters  by  Equation 
(2-1) .   The  next  step  in  the  calculation  is  to  apply  the  correc- 
tions to  the  intrusion  data.   The  corrections  required  during  low 
pressure  intrusion  are  to  account  for  the  fact  that  part  of  the  intru- 
sion recorded  is  actually  the  compression  of  the  air  that  was  trapped 
in  the  penetrometer  during  the  mercury-filling  period  (refer  to  Appen- 
dix A-2) .   The  high  pressure  correction  to  the  intrusion  data  is  to 
account  for  the  compressibility  of  the  sample,  penetrometer  and  mercury, 
as  well  as  the  temperature  changes  that  occur  as  a  result  of  pressuri- 
zation.   A  "blank  test"  for  this  correction  had  been  outlined  in  a  pro- 
posed ASTM  standard,  subcommittee  D18.14.   In  this  study,  a  correction 
of  0.001  ml  per  ml  of  mercury  in  the  penetrometer  at  103,500  kPa 
(15,000  psi),  as  noted  by  Bhasin  (1975),  was  applied. 

A  computer  program  developed  by  Reed  (1977)  was  utilized  to 
reduce  the  PSD  data.   This  program  consisted  of  the  above  calculations 
and  corrections.   A  significant  modification  of  this  program  has 
been  made,  especially  for  the  graphical  presentations  of  the  PSD 
data  and  the  finite  difference  approximation  (see  Section  4.3.2).   A 
listing  of  this  modified  computer  program  and  sample  output  are  given 


52 


in  Appendix  A-3. 

4.3.2   Data  Presentation 

The  PSD  data  are  normally  expressed  in  terms  of  the  cumulative 
volume  of  pore  space  intruded  as  a  function  of  pore  diameter.   Since 
a  wide  range  (e.g.,  several  orders  of  magnitude)  of  pore  sizes  is 
measured,  it  is  often  convenient  to  present  pore  diameter  on  a  log- 
arithmic scale.   A  graphical  cumulative  distribution  with  the  intruded 
volumes  per  gram  on  the  ordinate  to  an  arithmetic  scale  and  the  limi- 
ting pore  diameters  (  called  apparent  pore  diameters  )  on  the  abscissa 
to  a  logarithmic  scale  is  a  common  way  to  present  the  PSD  test  results  , 

e.g..  Diamond  (1970),  Ahmed  (1971),  Bhasin  (1975),  Cebeci,  et  al. 
(1978),  and   the   proposed  ASTM  standard,  subcommittee  D18.14.   From 
Equation  (2-5a)  it  is  evident  that  the  ordinate  of  the  density  function 
plot  is  -d(V  -  V)/dx.  Multiplying  this  derivative  quantity  by  x/2.303 
and  rearranging  gives  -d(V  -  V)/d(log  x).  Since  the  cumulative  dis- 
tribution is  expressed  as  (V  -  V)  versus  log  x,  it  is  preferred  to  plot 
the  quantity  -d(V  -  V)/d(log  x)  instead  of  -d(V   -  V) /dx  versus  log  x. 
The  density  function  value  can  be  obtained  by  (1)  directly  taking  the 

derivative  of  F    (x)  which  is  a  curve-fitted  function  of  the  above 
vx 

cumulative  pore  size  distribution  curve,  or  (2)  using  a  finite  differ- 
ence approximation.  The  first  method  is  obvious  and  its  accuracy  de- 
pends on  goodness  of  curve-fitting.  In  this  study,  the  pore  size  den- 
sity function  of  soils  is  obtained  in  this  manner  (see  Section  4.3.3). 
The  second  method  was  first  used  by  Garcia-Bengochea  (19  78)  for  the 
graphical  presentation  of  the  PSD  test  data.  However,  since  pore  size 
Is  a  continuous  random  variable,  it  is  appropriate  to  present  a  pore 


size  density  function  instead  of  a  frequency  polygon,  which  was  used 
by  Garcia-Bengochea .   To  approximate  the  density  function  a  frequency 
histogram  is  plotted  (i.e.,  A(V  -  V)  versus  log  x,  see  Figure  4-3a) . 
For  a  good  approximation,  the  class  width  ,  A(log  x) ,  has  to  be  an 
appropriate  constant. 

Assume  a  constant  c  is  chosen,  then 

log  x^_^  -  log  x^  =  c  .  (4-5) 

or  -^^-^  =  10  (4-6) 

^i 

From  Equation  (2-1) ,  it  is  evident  that  the  applied  pressure  is 
inversely  proportional  to  the  pore  diameter.   Therefore  this  approxi- 
mation method  requires  that  the  applied  pressure,  P  be  raised  as  follows: 

P. 
1       c 

=  10  (4-7) 


^-1 

Note  that  A(V  -  V)  represents  change  in  intruded  pore  volume, 
which  is  the  difference  between  the  two  successive  intrusions,  and  it 
can  be  calculated  directly  from  intrusion  record.   The  pore  size  den- 
sity function  plot  is  approximated  by  connecting  the  points  of  Ll^g 
d.,  f(log  d.)J  as  shown  in  Figure  4-3b ,  where  log  d.,  the  midpoint  of 
each  class,  is  calculated  as  follows: 

log  d^  =  I  (log  x._^  +  log  x^)   "  (4-8) 

or  d  =  /  X    •  X.  (4-9) 

and  f  (log  d.)  is  the  probability  density  at  log  d.,  which  is  calcu- 
lated as  follows: 
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A  (V   -  V) 

f(log  d^)  =  ^ 

A  (log  x) 


(V^  -  V)  at  P.  -  (V^  -  V)  at  P._^    ^,_^q) 


t*   "v 


or 


log  X.    -  lo^ 
fdog  d^)  .  ^(v  -  V)/c  .     ■         (4-11) 


A  few  specimens  were  tested  and  presented  according  to  this 
finite  difference  approximation  method.   Figure  4-4  shows  a  compari- 
son of  the  above  two  approaches.   It  is  evident  that  these  two  methods 
match  up  very  well  if  an  appropriate  constant   c  is  used.   The  value 
of  c  should  be  as  small  as  possible.   However  the  appropriate  value  -^f 
c  depends  on  the  pore  size  characteristics  of  soils.   The  wider  the 
pore  size  range  is  the  larger  Che  value  of  c  should  be.   In  this  stuay, 
the  value  of  c  varied  from  0.1  to  0.2. 

In  addition,  some  investigators  (e.g.,  Lohnes ,  et  al.,  197f.  and 
Diamond,  et  al . ,  1972)  prefer  to  express  pore  size  distribution  as  a 
dimensionless  term  (e.g.,  percentage)  for  the  convenience  of  che  mathe- 
matical treatment.   Equations  (2-6)  and  (2-7)  provide  the  theoretical 
basis  for  this  approach.   If  the  pore  size  distribution  data  are  to  be 
correlated  with  soil  properties,  it  is  desirable  to  express  the  pore 
size  distribution  in  mathematical  forms.   Hence  it  is  convenient  to 
express  pore  size  distribution  as  a  percentage.   In  Equation  (2-7), 
the  term  V   represents  the  total  void  volume  that  can  be  intruded  by 
mercury.   According  to  this  definition,  the  cumulative  distribution 
for  pores  with  size  greater  than  a  lower  limit  pore  diameter  (which 
depends  on  the  capacity  of  porosimeter)  should  be  100%  no  matter  tiow 
much  unintruded  pore  space  is  left.   On  the  other  hand,  if  the  term  V 
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The  unit  of  density  function  is  indeed  (cc/g/chosen  diameter  interval), 
but  the  chosen  diameter  interval  is  dimensionless  (see  Eq .  4-11). 
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represents  total  void  volume  that  exists  in  the  sample,  the  cumulative 
distribution  for  pores  with  size  greater  than  a  lower  limit  pore  dia- 
meter may  not  be  100%,  but  depends  on  how  much  pore  space  is  unintruded. 

4.3.3  Pore  Size  Distribution  Curves  and  Fittings  of  the  Curves 
Following  the  preceeding  discussion  of  data  presentation,  the  pore 
size  distribution  curves  directly  generated  from  the  PSD  data  are 
shown  in  Appendix  A-A.   In  addition  to  the  cumulative  distribution 
curves,  some  figures  also  show  the  density  function  curves  whenever 
the  finite  difference  approximation  method  is  applicable. 

At  least  three  PSD  specimens  were  tested  for  each  compacted  sample 
and  two  compacted  samples  were  tested  for  each  sample  code.   In  general, 
the  replication  of  the  PSD  tests  was  satisfactory.   However  some  scatter 
is  present  in  these  data  points.   The  reasons  for  scatter  are: 

1.  Some  of  the  specimens  might  have  suffered  more  disturbance 
than  others  during  trimming  and  transportation  of  the  speci- 
mens.  This  is  especially  a  problem  for  Soil  B  in  which  in- 
sufficient cohesion  is  created  by  the  small  amount  of  clay 
(i.e.,  10%)  . 

2.  The  compacted  samples  themselves  might  not  be  completely 
homogeneous . 

3.  Measurement  errors  may  also  account  for  scatter  in  replicate 
pore  size  distribution  curves,  especially  when  a  trapped  air 
pocket  on  the  surface  of  the  specimen  happens  at  the  filling 
stage  (see  Appendix  A-2)  or  when  the  pore  volume  of  a  specimen 
exceeds  the  intrusion  capacity  of  the  porosimeter. 

4c  The  variability  of  replicate  compacted  samples  may  also  account 
for  the  scatter  shown  in  these  figures. 


58 


A  plecewise  cubic  function  based  on  the  least  squares  method  was 
selected  to  fit  the  PSD  curve.   A  coii.iutor  program  called  "SPLINE"  in  the 
Purdue  University  CDC  6500-6600  computer  library  was  used  to  do 
the  least  squares  analysis.   A  description  of  "SPLINU"  and  some  sample 
output  are  listed  in  Appendix  A-5.   The  coefficients  of  the  piecev/ise 
cubic  functions  for  all  samples  are  also  listed  in  this  Appendix. 

The  next  step  is  to  calculate  the  functional  values  of  the  pore 
size  distribution  and  its  probability  density  and  then  to  present  the 
graphical  results  of  the  PSD  tei.ts.   Figure  4-5  shows  a  flow  chart  of 
the  procedures.   A  computer  program  for  this  calculation  is  listed  in 
Appendix  A-6.   Figure   4-6  through  Figure  4-15  show  these  graphical 
results.   Each  of  the  figures  consists  of  a  cumulative  distribution 
and  its  density  function   with  the  volume  for  gram  on  the 
ordinate  to  an  arithmetic  scale  and  the  limiting  pore  diameter  (or 
called  apparent  pore  diameter)  on  the  abscissa  to  a  logrithmic  scale. 

The  density  function  provides  a  useful  complement  to  the  cumula- 
tive distribution  curve  for  representing  pore  size  characteristics. 
The  pore  size  density  function  curves  displayed  bimodal  characteristics 
for  clay-sand  mixtures  (i.e..  Soils  B  and  C) .   This  was  the  case  found 
by  Garcia-Bengochea  (1978)  for  clay-silt  mixtures.   However,  a  single 
mode  characteristic  was  found  for  Soil  A.  This  is  logical  since  Soil 
A  contains  only  sand  grains. 

For  the  bimodal  characteristics  of  Soils  B  and  C,  the  small  pore 
mode  occurred  consistently  around  0.1  ym  and  the  large  pore  mode  was 
between  10  and  100  ym.   For  the  single  modal  characteristic  of  Soil 
A,  the  large  pore  mode  occurred  around  200  pm.   It  is  obvious  that  the 
mode  and  its  probability  density  may  serve  as  useful  qualitative 
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' 

' 

Calculation  of  Functional 
Values  from  SPLINE 
(Appendix  A-6  ) 

^ 

PSD  Curves 

Figure  4-5 


Flow  Chart  of  Procedures  for  the 
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Figure  4-6     PSD  Curves  of  SCIMW 
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Figure  4-7    PSD  Curves  of  SCIMO 
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Figure 4-8   PSD  Curves  of  SCIMD 
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Figure  4-9     PSD  Curves  of  SC3MW 
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Figure  4-10     PSD  Curves  of  SC3M0 
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Figure  4-11    PSD  Curves  of  SC3MD 
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Figure  4-12     PSD  Curves  of  SC3H0-D 
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Figure  4-13    PSD  Curves  of  SC3HD 
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Figure  4-14     PSD  Curves  of  OPL 
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Figure  4-15     PSD  Curves  of  OPH 
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descriptors  for  the  soil  fabrico 

For  the  purpose  of  correlating  the  theoretical  permeability 
model  with  the  pore  size  distribution  data,  it  is  preferable  to  ex- 
press the  cumulative  distribution  and  density  function  of  pore  size 
as  a  percentage.   It  should  be  noted  that  at  the  maximum  pressure 
exerted,  about  10  to  15  percent  of  the  pore  space  was  not  intruded 
for  the  soils  tested.   The  reasons  for  this  have  not  been  established 
with  certainty.   Some  of  the  pores  may  have  entry  ways  too  narrow  to 
penetrate  with  mercury  at  103,500  kPa  (15,000  psi) ,  and  most  likely 
some  of  them  may  be  in  a  network  of  pores  not  continuous  with  the 
boundaries  of  the  specimen.   Moreover,  a  minor  shrinkage  of  the 
specimens  during  sample  drying  could  account  for  part  of  the 
discrepancy. 

A  significant  unintruded  pore  space  is  not  unusual  for  pore  size  dis- 
tributions determined  by  the  mercury  intrusion  technique.   About  25%  of 
the  pore  space  was  not  intruded  by  mercury  in  the  study  by  Diamond,  et  al. 
(1972)  for  hydrated  cement  paste.   A  value  of  20%  unintruded  pore  space 
was  reported  by  Lohnes,et  al.  (1976)  for  selected  Hawaiian  soils  (e.g., 
Manana  soil,  consisting  of  70%  sand,  16%  silt  and  14%  clay).  A  45%  un- 
intruded pore  space  was  reported  by  White  (  1980)  for  a  medium  plastic 
clay,  St.  Croix  clay.   However,  only  about  5%  unintruded  pore  space 
was  reported  by  Garcia-Bengochea  (1978)  and  Reed  (1977)  for  clay-silt 
mixtures.   Note  that  the  above  values  are  all  under  the  maximum  intruded 
pressure  of  103,500  kPa  (15,000  psi)  except  the  one  used  by  Diamond, et  al. 
(1972),  in  which  a  maximum  pressure  of  517,500  kPa  (50,000  psi)  had 
been  applied.   With  these  facts  in  mind,  the  discontinuity  of  pores 
may  be  the  main  reason  of  the  existence  of  the  unintruded  pore  space 
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if  a  proper  specimen  drying  method  is  used.   However,  Reed,  et  al. 

( 1980)  have  expressed  the  view  that  nearly  all  pore  space  in  soil  systems 

is  connected. 

Fortunately,  this  unintruded  pore  space  effect  is  of  minor  impor- 
tance in  this  study  which  deals  with  flow  through  porous  materials,  be- 
cause the  pore  channels  through  which  water  passes  are  probably  larger 
than  those  penetrated  by  mercury  at  pressure  of  15,000  ps i . 

An  attempt  at  a  correction  for  this  unintruded  pore  space  effect 
was  made  by  Meyer  (1953).   A  modified  form  of  his  work  is  given  in 
Appendix  A-7 

Because  of  this  unintruded  pore  space,  the  graphical  results  of 
the  pore  size  distribution  and  its  density  function  are  presented  in 
both  manners  discussed  in  Section  4.3.2.   An  example  is  sho^^m  in 
Figures  4-16  and  4-17.   The  rest  of  them  are  shown  in  Appendix  A-8 . 

4.3.4   Comparison  of  Pore  Size  Curves 

Figures  4-18,  4-19,  and  4-20  show  the  effect  of  water  content  on 
pore  size  distribution.   For  SC3M  samples  (see  Figure  4-18)  ,  the  pore 
size  density  function  shows  a  bimodal  characteristic  with  differences 
in  the  large  pore  mode,  and  almost  no  changes  in  the  small  pore  mode. 
The  large  pore  mode  decreases  as  water  content  increases  from  dry  of 
optimum  to  optimum  water  content.   There  is  little  difference  in  the 
large  pore  mode  of  the  optimum  and  the  wet  side  samples.   In  general, 
the  trend  as  indicated  above  is  similar  to  that  found  by  Garcia- 
Bengochea  (1978) . 

Comparing  the  small  pore  modes  displayed  in  Figures  4-18  and  4-19, 
it  is  evident  that  the  small  pore  mode  is  almost  the  same  in  spite  of 
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Figure  4-16     PSD  Curves  of  SC3MD  ( Percent  of 
Total  Pore  Volume  Existing) 


73 


100 


UJ 

n_ 


CO 


30 


60 


40 


20 


\ 

\ 

\ 

' 

D 

^ 

^^ 

r-^--^. 

N 

\ 

\ 

W„ 

.01  0.1  1.0  10.0  100.0  1000 

LIMITING  PORE  DIfl.METER    IN  MICRONS 


U 

UJ 


UJ 

a 


HUU 

LESEND 

■      0     SflKPLE  CODE 

SC3« 

3nn 

- 

?nn 

inn 

/ 

\ 

A 

0 

V 

,.,^ 

y  ^ 

L 

.01     0.1      1.0      10.0     100.0 
LIMITING  PORE  DIRMETER  IN  MICRONS 


1000 


Figure  4-17     PSD  Curves  of  SC3MD  (Percent  of 
Total  Pore  Volume  Intruded  ) 
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Figure  4-18    Comparison  of  PSD  Curves  for  Samples 
Compacted  at  Different  Water  Contents 
(SC3M  series) 
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Figure  4-19    Comparison  of  PSD  Curves  for  Samples 
Compacted  at  Different  Water  Contents 
(SC3H  series) 
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Figure  4-20    Comparison  of  PSD  Curves  for  Samples 
Compacted  at  Different  Water  Contents 
(SCIM  series) 
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an  increase  of  the  compactive  effort.   This  phenomenon  will  be  dis- 
cussed later.   Unlike  the  trend  in  Figure  A-18,  the  large  pore  modes 
of  the  dry  side  and  the  optimum  samples  are  quite  similar  in  Figure 
4-19,  although  they  do  show  a  more  significant  difference  in  the  cu- 
mulative distribution  curves.   This  may  be  due  to  inaccurate  location 
of  points  of  the  optimum  samples  (see  Figure  4-1) .   The  optimum 
samples  were  indeed  slightly  dry  of  optimum  water  content,  thus  the 
samples  were  coded  as  SC3H0-D.   Comparing  results  of  permeability 
trsts  (see  Table  4-6) ,  the  difference  between  permeability  of  Sample 
SC3HD  and  Sample  SC3H0-D  is  much  smaller  than  that  of  Sample  SC3MD 
and  Sample  SC3M0.   The  similar  trend  when  comparing  pore  size  charac- 
teristics and  those  of  permeability  implies  that  pore  size  distri- 
bution could  be  a  sensitive  indicator  of  permeability. 

Figure  4-20  shows  that  the  difference  between  por^  size  distribut- 
ion of  the  dry  side  and  optimum  or  wet  side  samples  is  small.   The 
reasons  noted  in  Section  4-2  for  comparison  of  permeability  measure- 
ments on  the  same  soil  compacted  at  different  water  contents  may  also 
contribute  to  the  small  differences  between  these  pore  size  character- 
istics. -   •      , 

As  noted  in  Section  2.1.3,  seepage  occurs  only  through  the  pores, 
and  what  is  directly  related  to  the  permeability  is  the  characteristics 
of  the  pores  themselves.   Thus  the  permeability  of  each  of  the  samples 
increases  with  increasing  large  pore  mode.   Insignificant  difference 
in  the  large  pore  mode  is  reflected  by  insignificant  difference   in 
permeability. 

Figure  4-21  and  Figure  4-22  show  the  effect  of  compactive  effort 
on  pore  size  characteristics.   Note  that  these  comparisons  are  only  semi- 
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Figure  4-22    Comparison  of  PSD  Curves  for  Samples 
Compacted  at  Different  Compactive 
Efforts  (Soil  A) 
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qualitative  because  the  water  content  was  not  controlled  at  the  same 
value.   Figure  4-21  shows  that  as  the  compactive  effort  was  increased, 
the  large  pore  mode  decreased,  while  the  small  pore  mode  remains  un- 
changed.  Figure  4-22  also  shows  the  same  trend  except  that  there  is 
no  small  pore  mode.   Recalling  the  permeability  of  these  soils,  again 
it  decreases  with  decreasing  large  pore  mode.   This  trend  was  also 
found  by  Garcia-Bengochea  (1978). 

Figures  4-23,  4-24  and  4-25  compare  the  pore  size  characteristics 
of  Soils  B  and  C  compacted  by  medium  effort  dry  of  optimum,  optimum 
and  wet  of  optimum  water  contents,  respectively.   The  large  pore  mode 
of  the  soils  decreases  with  increasing  clay  content,  while  the  small 
pore  mode  remains  constant,  but  its  probability  density  increases  with 
increasing  clay  content. 

Figure  4-26  compares  the  pore  size  characteristics  of  Soils  A, 
B,  and  C  with  different  compaction  variables.   The  small  pore  mode 
remains  the  same  in  spite  of  differences  in  compactive  effort  and 
amount  of  clay,  while  its  probability  density  decreases  with  decreas- 
ing clay  content.   Note  that  there  is  no  small  pore  mode  in  soil  A, 
which  follows  from  the  fact  that  Soil  A  has  no  clay  content.   The 
large  pore  mode  of  the  soils  increases  with  increasing  sand  content 
as  mentioned  above.   Comparing  the  permeability  values  and  pore  size 
characteristics  of  these  soils  (see  Table  4-6  and  Figure  4-26)  shows 
that  the  permeability  decreases  with  a  decrease  in  the  large  pore  mode. 

For  the  soils  tested,  the  large  pore  mode  and  its  probability 
density  appears  to  be  the  best  indicator  of  changes  in  soil  fabric 
caused  by  varying  water  content,  compactive  effort,  and  soil  type. 
This  was  also  found  by  Garcia-Bengochea  (1978) . 
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Figure  4-23    Comparison  of  PSD  Curves  for  Different 
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Figure  4-24    Comparison  of  PSD  Curves  for  Different 
Soils  (B and  C )  Connpacted  by  Medium 
Effort  at  Optimum  Water  Content 
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Figure  4-25    Comparison  of  PSD  Curves  for  Different 
Soils  (B  and  C)  Compacted  by  Medium 
Effort  at  Wet  of  Optimum  Water  Content 


84 


u> 


o 

o 


t— • 

=3 
BC 

t— 


.01     0.1      1.0      10.0      100.0 
LIMITING  PORE  DIAMETER  IN  MICRONS 


1000 


1.0 


CO 


Q    .SO 

llJ 
o 


t,  .25 

UJ 

o 


cr 
o 


.00 


LEGEND 

•  0     SAMPLE  CCDE       SC3M0 

•  A     SflffLE  CODE       SC3HD 

--     +     SHMPLE  CODE       SCIMD 
■  ■      X     Sfr?L£  CODE       OPL 
•'     ♦     SflMPLE  CCDE       OPH 

h 

<  • 

< 
: 

i 

•      .Ji^ 

Va=^«^ 

j>A 

is. 

.01  0.1  1.0  10.0  100.0 

LIMITING  PORE  DIAMETER    IN  MICRONS 


1000 


Figure  4-26    Comparison  of  PSD  Curves  for  Different 
Soils  ( A, B  and C) 


85 


The  hypothesis,  proposed  by  Garcia-Bengochea  (19  78),  that  the 
shape  of  the  differential  pore  size  distribution  (called  pore  size 
density  function  in  present  study)  in  the  small  pore  region  (referred 
to  as  small  pore  mode)  may  serve  as  a  "signature"  of  a  given  clay  type, 
is  supported  by  the  data  in  this  study.   Again,  the  small  pore  mode 
is  related  to  the  size  range  of  clay  particles  present,  and  the  pro- 
bability density  of  the  small  pore  mode  is  in  proportion  to  the  amount 
of  clay  minerals  in  the  soil. 

4.4  PREDICTION  OF  PERMEABILITY  BY  PROPOSED  MODEL 

The  objective  of  this  section  is  to  correlate  the  permeability  and 
the  PSD  of  the  selected  soils  by  applying  the  theoretical  model  pro- 
posed in  Section  2.2.2. 

Before  going  into  a  detailed  discussion  of  the  proposed  model,  it 
is  worthwhile  to  examine  the  following  functional  relationships  between 
void  ratio,  e  and  permeability,  k  as  presented  by  Taylor  (1948)  and 

Lambe  (1951) : 

3 

(a)  k  vs .  e  /(I  +   e) 

(b)  k  vs.  e^/(l  +  e) 

(c)  k  vs .  e 

(d)  log  k  vs .  e 

Since  the  permeabilities  of  the  soils  studied  vary  over  five  orders  of 
magnitude,  it  is  very  unlikely  that  linear  relationships  for  these 
relationships  exist.   Therefore  the  logarithm  of  k,  instead  of  k,  was 
used  for  the  regression  analyses.   The  statistical  programs  of  the 
Statistical  Package  for  the  Social  Sciences  (SPSS)  in  the  Purdue  Univer- 
sity CDC  6500-6600  computer  library  were  used  to  do  least  squares  ana- 
lysis and  regressions  in  this  study.   The  coefficient  of  determination 
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2 
(R  )  of  the  regressions  of  k.  with  respect  to  the  functional  parameters 

2 
of  e  mentioned  above  are  shown  in  Table  4-7.   Although  R  was  around 

0.6,  the  differences  between  the  measured  and  the  calculated  permea- 
bilities are  generally  greater  than  one  order  of  magnitude  and,  in  some 
samples  are  greater  than  two  orders  of  magnitude.   The  results  noted 
above  imply  that  gross  volumetric  parameters  (e.g.  void  ratio  or  poro- 
sity and  their  functions)  cannot  appropriately  account  for  changes  in 
permeability  of  compacted  soils. 

Previous  discussions  in  Section  2.2  provided  the  theoretical 
background  for  predicting  the  permeability  from  the  PSD  of  soils .   From 
the  qualitative  observations  discussed  in  Section  4.3.4,  it  is  evident 
that  the  permeability  is  related  to  the  distribution  of  pore  sizes  about 
the  large  pore  mode.   Moreover,  the  proposed  model  is  very  sensitive  to 
the  larger  pore  diameters  (see  Equation  2-19) .   Therefore  the  pores  of 
diameter  less  than  lym  were  considered  to  be  negligible  in  the  calcu- 
lation of  permeability.   A  physical  explanation  of  choosing  lym  as  the 
lower  limit  of  pore  sizes  which  is  accessible  to  flow  is  the  concept 
of  "threshold  diameter"  given  by  Murota  and  Sato  (1969).   In  their 
study,  a  threshold  pore  diameter  of  1.8ym  was  chosen  for  silty  clay. 
On  the  other  hand,  the  upper  limit  of  pore  sizes  in  this  study  was 
bounded  by  the  pressure  capacity  of  the  mercury  porosimetry.   Hence 
Equation  (2-19)  can  be  rewritten  as: 

2  f   ^max  f  ^max 
k  =  ^         /       X  G(x.)  .  f(x.)  •  f(x.)  dx.  dx.       (4-12) 
^  Jl      7l  ^        ^       J     1    J 

Numerical  integration  of  Equation  (4-12)  was  accomplished  by  a 

computer  program  listed  in  Appendix  A-9 .   Note  that  the  equivalent 

diameter,  x  was  assumed  to  be  the  smaller  of  each  pair  of  connected 
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pores.   The  term  G(x.)  defined  in  Equation  (2-18)  contains  a  undeter- 
mined function,  V.   Thus  in  order  to  predict  permeability  from  Equation 
(4-12) 5  V  must  be  assumed.   In  this  study,  V  was  assumed  to  be  a 
linear  function  of  pore  diameter  (i.e.  V(x.)  =  v  •  x.)  in  order  to 
simplify  the  calculations.  The  pore  size  density  function  f  was  deter- 
mined previously  (see  Section  4.3.3  and  Appendix  A-5)  and  was  expressed 
as  a  dimensionless  term.   In  other  words,  f  was  a  piecewise  cubic  func- 
tion which  fitted  its  corresponding  PSD  curve  (i.e.,  Figures  A-8-1 
through  A-8-20) . 

The  computer  program  listed  in  Appendix  A-9  provides  an  iterative 
method  to  seek  a  suitable  value  of  v  at  whuch  a  good  prediction  of 
permeability  is  obtained.    In  this  computer  program  the  criterion  of 
a  good  prediction  of  permeability  is  determined  by  the  user.   In  the 
present  study,  a  criterion  of  a  15%  error  with  respect  to  the  measured 
permeability  was  arbitarily  chosen.   Two  other  criteria  were  number 
of  iterations  and  the  relative  change  of  the  calculated  permeability 
between  the  two  successive  values  of  v.   For  each  iterative  run,  the 
main  feature  of  the  computer  program  listed  is  the  numerical  integra- 
tion of  the  double  integral  in  Equation  (4-12) .   A  subroutine  called 
"CADRE"  in  the  Purdue  University  CDC  6500-6600  computer  library  was 
utilized  to  carry  out  this  numerical  integration. 

The  output  of  the  interations  discussed  above  for  the  samples 
tested  is  also  given  in  Appendix  A-9. 

Following  the  above  discussions  and  the  results  given  in  Appendix 
A-9,  the  values  of  v  for  all  samples  studied  are  listed  in  Table 
4-8.   The  next  step  is  intended  to  correlate  v  to  the  pore  size 
characteristics  of  soils.   Table  4-9  lists  some  of  the  possible 


TABLE  4-8   The  v  -  Values 
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Sample 
Code 


from  iteration 


from  Equation  (A-13) 


OPL 

OPH 
SC3MW 
SC3M0 
SC3MD 
SCIMW 
SCIMO 

SCIMD 
SC3H0-D 

SC3HD 


0.1500 
0.1500 
0.0020 
0.0012 
0.0350 
0.0600 
0.0700 
0.4000 
0.0025 
0.0070 


0.1400 
0.1360 
0.0014 
0.0015 
0.0386 
0.0510 
0.1500 
0.2420 
0.0027 
0.0069 


TABLE  4-9   Possible  Correlations  between  PSD 
Parameters  and  v 
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Dependent  Variable 


Independent  Variable 
(PSD  parameters) 


V 

or 

ln(v) 

V 

or 

ln(v) 

V 

or 

ln(v) 

V 

or 

ln(v) 

V 

or 

ln(v) 

V 

or 

ln(v) 

V 

or 

ln(v) 

V 

or 

ln(v) 

V 

or 

ln(v) 

V 

or 

ln(v) 

V 

or 

ln(v) 

V 

or 

ln(v) 

V 

or 

ln(v) 

V 

or 

ln(v) 

V 

or 

ln(v) 

''25  (see  remarks) 

75 

^^^25  (see  remarks) 
LD, 


D 


LD 


50 

75 


•^25'  °50'  °75 

LD23.  ^^0'  ^°75 

^0''°80 

io/(°10  ^  °60^ 

°60/°10 

^0/°25 

^40/°10 

°60/°10'  ^o/°25' 

°60/°25'  ^0^%' 


°40/°10 
^25/^10 


Remarks 


a)   D  3  means  25  Percentile  (a  diameter)  of  the  cumula- 
tive distribution  curve,  i.e.,  the  subscripts  re- 
present the  percentile  of  the  distribution. 


b)   LD   means  the  log  of  diameter  D  ^ 
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correlations  between  PSD  parameters  and  v-values.   Simple  linear  re- 
gression by  the  method  of  least  squares  was  used  to  evaluate  these 

possible  correlations,  and  the  following  satisfactory  linear  relation 

2  2 

was  established  with  R  =  0.97  and  an  adjusted  R  =  0.96: 

■^       a 

In  V  =  -  6.15  -  40.07  (-^)  +  118.79  (-^)  -  78.35  (-^)       (4-13) 

°25  10  10 

With  Equation  (A-13)  and  the  PSD  of  soils,  the  value  of  v  can  be 
determined  (see  Table  4-8).   The  permeability  can  then  be  calculated 

b>  using  the  computer  program  described  above.   Table  4-10  compares 
the  calculated  and  the  measured  permeabilities.   A  linear  regression 
was  performed  and  the  result  as  shown  in  Figure  4-27  was  excellent. 
In  summary,  the  proposed  permeability  model  can  provide  an  ex- 
cellent predictive  tool.   Thus  a  functional  relationship  between  the 
PSD  and  the  permeability  of  compacted  soils  is  clearly  demonstrated. 

4.5   APPLICABILITY  OF  THE  PROPOSED  PERMEABILITY  MODEL 

Attempts  were  made  to  broaden  the  range  of  prediction  of  per- 
meability using  independently  determined  data.   Ten  sets  of  PSD 
curves  along  with  permeability  values  were  taken  from  Garcia- 
Bengochea  (1978) .   The  procedures  described  above  were  carried  out 

for  both  the  author's  and  Garcia-Bengochea 's  data.   The  following 

2 
regression  equation  was  established  for  the  v-value  (R  =  0.89): 

In  V  =  _  3.6031  -  2.5957  (-^)  +  14.9678  (—-) 

25  10 

-  7.2562  (^)  +  4.6320  (-^)  -  2.7129  (--^) 

10  "25  "lO 

-  3.8473  (-^)  (4-14) 

'^40 


92 


TABLE  4-10   Comparison  of  Predicted  and  Measured  Permeability 


Sample  Code 


Predicted  Permeability 
(cm/sec) 


Measured  Permeability 
(cm/sec) 


OPL 

OPH 
SC3MW 
SC3M0 
SC3MD 
SCIMW 
SCIMO 
SCIMD 
SC3H0-D 
SC3HD 


2,44  X  10 


-2 


1.43  X  10 


4.40  X  10 


-2 


-7 


9.85  X  10 


4.85  X  10 


4.98  X  10 


2.24  X  10 


-4 


-3 


2.40  X  10 


-3 


7.31  X  10 


1.70  X  10 


-6 


2.5  X  10 


1.4  X  10 


-2 


-2 


6.0  X  10 


-7 


7.3  X  10 


-7 


4.5  X  10 


-5 


5.3  X  10 


-4 


1.3  X  10 


-3 


3.5  X  10 


-3 


6.1  X  10 


1.9  X  10 


-7 
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Figure  4-27    Prediction  of  Permeability  by  Proposed  Model 
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The  values  of  v  calculated  from  Equation  (4-lA)  were  then  used 
In  the  calculation  of  permeability.   The  information  from  each  data 
set,  including  sample  code,  PSD  curve-fitting  result,  the  v-value, 
and  the  calculated  and  measured  permeabilities  are  given  in  Appendix 
A- 10.   Note  that  the  computer  program  listed  in  Appendix  A-9  was  used 
to  perform  the  above  calculations.   In  addition,  a  regression  analysis 
was  performed  on  the  permeability  data  (i.e.,  the  calculated  and  the 
measured  permeability).   The  result  is  shown  in  Figure  4-28.   It  is 
evident  that  the  proposed  permeability  model  provides  as  excellent 
theoretical  tool  for  the  prediction  of  permeability  for  a  wide  range 
of  soil  types. 
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Regression  Equation  (Rg  ■  .98) 

Y-  .99  X-.03 
where 

Y"log  of  permeabilily (calc. 
X»  log  of  permeability (meas.) 


10               10 
Permeability  (cm/sec) Measured 


Figure  4-28    Prediction  of  Permeability  by  Proposed 
Model  (Author's  and  Garcia-Bengochea's 
Data) 
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CHAPTER  V   SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 

5.1   SUMMARY 

The  theoretical  background  for  the  determination  of  the  pore  size 
distribution  of  soils  has  been  critically  reviewed.   Mathematical 
interpretations  of  PSD  were  given  in  detail.   For  the  soils  studied, 
both  the  pore  size  distributions  and  pore  size  density  functions  were 
presented.   The  pore  size  distributions  from  original  PSD  data  were 
expressed  as  intruded  volume  per  gram  of  specimen.   In  addition,  they 
were  also  expressed  as  a  dimensionless  term  (e.g.,  percentage)  for 
convenient  mathematical  treatment.   The  pore  size  density  function 
could  be  obtained  either  from  direct  differentiation  of  a  curve-fitted 
cumulative  distribution  function,  or  by  finite  difference  approxima- 
tion.  Moreover,  the  results  of  both  methods  agree  very  well. 

The  soils  used  in  the  study  were  a  sand  and  two  sand-clay  mixtures, 
Pluvial  compaction  was  used  to  prepare  compacted  samples  of  sand  while 
kneading  compaction  was  used  to  prepare  samples  of  sand-clay  mixtures. 
The  testing  program  consisted  of  the  permeability  test  and  the  PSD 
test.   The  experimental  procedures  of  these  tests  for  sand-clay  mix- 
tures (called  Soils  B  and  C  in  this  study)  generally  followed  Garcia- 
Bengochea  (1978).   For  the  sand  (  called  Soil  A  ),  a  small  amount,  0.5% 
by  weight,  of  powdered  phenolic  resin  was  mixed  with  sand  before  com- 
paction.  Oven  drying  at  a  temperature  of  170  C  for  1   h  was  carried 
out  after  compaction  in  order  to  create  a  "cohesion.   This  technique 
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enabled  PSD  specimens  to  be  prepared  and  PSD  tests  to  be  conducted. 
Thus,  the  characterization  of  sand  fabrics  by  their  PSD  was  possible. 

Qualitative  descriptions  of  the  pore  size  characteristics  of  the 
soils  studied  were  given  through  their  pore  size  density  functions. 
When  displayed  on  a  log  diameter  scale,  the  pore  size  density  function 
of  Soil  A  showed  a  single  modal  characteristic  while  those  of  soils 
B  and  C  showed  bimodal  characteristics.   The  location  of  the  mode  (at 
which  a  peak  was  observed)  and  its  probability  density  could  serve  as 
a  useful  qualitative  descriptor.   For  the  bimodal  characteristics  of 
Soils  B  and  C,  the  small  pore  mode  occurred  consistently  around  0.1 
^m  and  the  large  pore  mode  was  between  10  and  100  um.   For  the 
single  modal  characteristic  of  Soil  A,  the  large  pore  mode  occurred 
around  200  ym. 

The  hypothesis,  proposed  by  Garcia-Bengochea  (1978),  that  the 
shape  of  the  differential  pore  size  distribution  (called  pore  size 
density  function  in  the  present  study)  in  small  pore  region  (referred 
to  as  small  pore  mode)  may  serve  as  a  "signature"  of  a  given  clay 
type,  is  supported  by  the  data  obtained  in  this  study.   In  addition, 
changes  in  soil  fabric  caused  by  varying  water  content  and  compac- 
tive  effort  were  generally  reflected  by  changes  in  large  pore  mode 
and  its  probability  density. 

The  relation  between  the  compaction  variables  and  the  permea- 
bility found  in  this  study  generally  agree  with  the  trends  found  by 
Lambe  (1954),  Mitchell,  et  al.  (1965)  and  Garcia-Bengochea  (1978). 
For  the  prediction  of  permeability,  the  relation  between  permeability 
and  pore  geometry  was  critically  reviewed.   A  new  permeability  model 
based  on  Child  and  Collis-George's  (1950)  model  was  proposed.   The 
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new  model  had  the  following  characteristics:   (1)  a  probabilistic 
basis,  (2)  no  arbitrary  matching  factor  required,  and  (3)  considera- 
tion of  a  tortuosity  factor.   The  prediction  of  permeability  by  the 

2 
proposed  model  was  excellent  (e.g.,  the  R   of  regression  of  predicted 

permeability  with  respect  to  measured  permeability  equaled  0.99). 

Data  obtained  independently  by  Garcia-Bengochea  (1978)  was  also  used 

to  examine  the  proposed  model.   The  excellent  result  Le.g.,  R"  =  0.98 

3. 

for  both  author's  and  Garcia-Bengochea 's  (1978)  data^  confirmed  the 
applicability  of  proposed  model. 

5.2   CONCLUSIONS 

lo   A  technique  for  the  preparation  of  PSD  specimens  of  cohesion- 
less  sands  has  been  developed. 

2.  Finite  difference  approximation  method  is  an  appropriate 
approach  for  determining  pore  size  density  functions.   The 
graphical  presentation  of  pore  size  density  function  obtained 
from  a  finite  difference  approximation  agrees  very  well  with 
that  obtained  from  direct  differentiation  of  the  curve-fitted 
cumulative  distribution  function. 

3.  The  curve  fitting  technique  using  the  computer  program  "SPLINE" 
for  the  PSD  data  of  soils  studied  was  apparently  satisfactory. 

4.  Both  pore  size  distribution  and  pore  size  density  function  are 
useful  for  characterizing  the  pore  size  characteristics  of  the 
soils  studied  in  this  research.   The  mode  and  its  probability 
density  of  pore  size  density  function  curves  served  <is 
successful  qualitative  descriptors,  especially  in  characteriz- 
ing the  changes  in  soil  fabric  caused  by  varying  compact- 
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ion  variables.   The  pore  size  parameters  (e.g.,  some  percen- 
tile and  their  functions)  obtained  from  pore  size  distribu- 
tions were  successfully  used  to  predict  the  constant  v  of  the 
proposed  permeability  model. 

5.  The  pore  size  density  function  of  Soil  A  (Ottawa  sand)  was 
found  to  be  singly  modal  on  a  log  diameter  scale.   The  mode 
occurred  around  200  ym. 

6.  The  pore  size  density  functions  of  Soils  B  and  C  (sand-clay 
mixtures)  were  found  to  be  bimodal  on  a  log  diameter  scale. 
A  large  pore  mode  occurred  between  25  ym  and  100  ym,  and 

a  small  pore  mode  occurred  around  0.1  ym. 
The  following  conclusions  (7  through  12)  verify  conclusions 
previously  reached  by  Garcia-Bengochea  (1978). 

7c   The  large  pore  mode  decreased  as  water  content  increased  from 

dry  of  optimum  to  optimum  water  content. 
8c   The  large  pore  mode  decreased  as  the  compactive  effort  in- 
creased. 
9.   The  large  pore  mode  of  the  soils  tested  decreased  with 
increasing  clay  content. 
10.   The  small  pore  mode  remained  the  same  in  spite  of  differences 
in  compactive  effort  and  amount  of  clay,  while  its  probabi- 
lity density  decreased  with  decreasing  clay  content, 
lie   For  the  soils  tested,  the  large  pore  mode  and  its  probability 
density  appear  to  be  the  best  indicators  of  changes  in  soil 
fabric  caused  by  varying  water  content,  compactive  effort, 
and  soil  type. 
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12.  The  hypothesis,  proposed  by  Garcia-Bengochea  (1978),  that 
the  shape  of  the  differential  pore  size  distribution  (called 
pore  size  density  function  in  the  present  study)  in  the  small 
pore  region  (referred  to  as  small  pore  mode)  may  serve  as  a 
"signature"  of  a  given  clay  type,  is  supported  by  this  study. 

The  following  conclusions  (13  through  18)  apply  to  the  proposed 

permeability  model. 

2 

13.  The  coefficient  of  determination  (R  )  of  the  conventional 

relation  between  permeability  k  and  void  ratio  e  (e  vs  log  k) 
was  0.54  for  the  soil  samples  tested.   This  result  showed  that 
k  can  not  be  accurately  predicted  from  total  porosity  or  functions 
thereof. 

14.  A  model  has  been  proposed  which  relates  PSD  to  soil  permeability 
which  is  as  follows: 


k  =  ^ 


2 

32y 


/    /     ^   •  G(x  )  •  f(x.)  •  f(x.)  .  dx.  dx. 
Jo      Jo  J        1        J       1    J 


X.   -  X.   „ 


where  G(x.)  =  e   ^   V 

and  V  =  V  .  X 

J 

15.   The  following  regression  equation  for  the  constant  v  was  ob- 

2 
tained  for  the  samples  tested  in  this  study  (R  =  0.97  and 

R^  =  0.96)  : 
a 

°40  °40 

In  V  =  -  6.15  -  40.07  (-^)  +  118.79  {-—-) 

°25  10 

-  78.35  (^p^) 
10 
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16.  The  excellent  predictive  result  by  the  proposed  model  is 

evident  from  the  following  regression  equation  for  the  samples 

2 

tested  in  this  study  (  R  =  0.99): 

a 

Y  =  0.997X  -  0.010 

when      Y  =  the  log  of  the  calculated  permeability 

X  =  the  log  of  the  measured  permeability 

17c  Independent  data  have  also  been  used  to  examine  the  proposed 

model.  The  following  regression  equation  for  the  constant  v 

was  obtained  for  the  author's  and  Garcia-Bengochea 's  (1978) 

data  (R^  =  0.89  and  R^  =  0,84): 

a 

In  v   =  -3.06  -  2.59  (tt^)  +  14.97  (-^) 

^25  10 

-7.26  (-^)  +  4.63  (-^)  -  2.71  (-^) 
"lO  25  10 

-3.85  M) 
40 

18.   The  applicability  of  the  proposed  model  is  evident  from  the 
following  regression  equation  for  the  author's  and  Garcia- 
Bengochea  's  (1978)  data  (R^  =  0.98): 

Y  -  0.99X  -  0.03 
5.3   RECOMMENDATIONS  FOR  FUTURE  RESEARCH 
Future  research  is  recommended  in  the  following  areas: 
1.  Prediction  of  permeability  by  the  proposed  model  should 
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be  extended  to:  (a)  other  compacted  conditions  of  soils  similiar 
to  the  ones  tested  here,  (h)    other  laboratory  compacted  soils 
and  (c)  field  compacted  soils. 

2.  An  investigation  of  the  sand  fabric  through  variation  in  PSD 
caused  by  different  compaction  methods  and/or  different  sands 

but  at  the  same  total  porosity  should  be  carried  out.   An  attempt 
could  then  be  made  to  relate  PSD  of  sands  to  their  liquefaction 
potential . 

3.  This  study  established  procedures  by  which  PSD  can  be  measured 

in  sand.   These  procedures  can  now  be  used  to  explain  and  predict 
shear  and  compressibility  benavior  of  sands. 

4.  Techniques  for  determining  PSD  in  coarse  sands  and  fine  gravels 
are  needed  to  exploit  this  parameter  in  the  prediction  of  behavior 
of  (a)  drainage  filters  and  (b)  highway  subbases. 
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